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Abstract 

The  phenomena  of  resonant  radiation  pressure  is 
studied  as  it  applies  to  the  reflection  of  two-level  atoms 
from  an  evanescent  wave.  The  basic  mechanism  is  examined  in 
detail.  Three  applications  are  proposed,  analyzed,  and 
compared  to  existing  methods  of  neutral  beam  control.  The 
first  application  is  a  beam  conducting  pipe  which  can  contain 
and  decelerate  a  beam  of  atoms.  The  second  application  is  a 
parabolic  nozzle  which  can  collimate  and  focus  a  beam  three 
orders  of  magnitude  better  than  comparable  resonant  focusing. 
The  third  application  is  a  macroscopic  trap  for  slow  atoms 
which  allows  laser-free  observation  of  '’lOOcm/s  atoms  for 
times  measured  in  seconds.  A  computer  code  which  tracks 
atomic  trajectories  through  multiple  reflections  is  included. 
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ELECTROMAGNETIC  MIRROR  DYNAMICS 


I.  Introduction 


Background 

Beams  of  neutral  particles  are  extremely  useful  for 
studying  many  physical  phenomena  such  as  natural  fluorescence 
and  neutral  particle  scattering.  Until  recently,  the  only 
effective  means  of  manipulating  neutral  beams  was  with  slits 
and  baffles,  or  (where  possible)  by  inducing  a  dipole  in  the 
atoms.  Both  of  these  methods  have  very  limited  applications. 

Resonant  radiation  pressure  has  been  used  in  the  last 
ten  years  for  limited  control  of  neutral  beams  (Refs  1-7). 
This  pressure  is  a  quantum-mechanical  force  which  occurs 
when  a  strong  laser  beam  strikes  the  beam  of  atoms.  The 
pressure  can  be  used  to  deflect  (Refs  1-3),  focus  (Ref  5), 
and  accelerate  (Refs  7-8)  a  beam  of  atoms.  Atomic  trapping 
schemes  have  been  proposed  (Refs  9-13)  but  not  demonstrated. 
All  of  these  applications  have  used  laser  beams  which  propa¬ 
gate  through  a  vacuum  and  directly  strike  the  atomic  beams. 

The  direct  contact  of  the  atomic  and  laser  beams 
imposes  two  limitations,  the  first  is  due  to  different  atomic 
speeds  in  the  beam,  which  causes  Doppler  shifts  and  greatly 
affects  the  strength  of  the  resonant  radiation  pressure. 


The  second  limitation  is  due  to  the  continuous  pressure 
needed  in  any  trapping  scheme.  This  means  that  the  atoms 
are  constantly  in  the  laser  beam  and  affected  by  it.  Obser¬ 
vation  of  trapped  atoms  unaffected  by  the  laser  is  highly 
desirable  for  very  accurate  spectroscopy  (Ref  7). 

A  new  arrangement  for  using  resonant  radiation 
pressure  has  been  developed  which  involves  the  total  internal 
reflection  of  the  laser  beam.  This  reflection  creates  a 
thin  ('•10  micron)  evanescent  wave  above  the  reflecting  sur¬ 
face  which  could  be  used  to  deflect  a  resonant  neutral  beam. 
The  atoms  would  only  be  in  the  radiation  for  a  short  time, 
and  this  would  allow  unaffected  observation  of  the  atoms  if 
a  proper  mirror  geometry  could  be  found.  The  mirror  could 
also  be  curved  and  arranged  to  focus  and  manipulate  a  beam. 

This  new  arrangement  has  many  possibilities,  but  this 
type  of  reflection  is  very  different  fronf  Standard  elastic 
reflection.  The  direction  of  the  reflected  atomic  beam 
depends  on  its  incident  speed  and  direction.  It  also 
depends  on  the  reflection  angle  of  the  laser  beam.  In  addi¬ 
tion,  the  strength  of  reflection  (how  much  of  the  beam  will 
reflect)  is  very  sensitive  to  the  speed  and  direction  of  the 
atomic  beam  and  the  frequency  and  direction  of  the  laser. 

No  systematic  investigation  of  these  factors  has  been  made. 


Problem  Statement 


Evanescent  resonant  reflection  is  investigated  in 
order  to  find  applications  in  the  control  of  neutral  atomic 
and  molecular  beams.  The  investigation  focuses  on  overcoming 
the  unusual  nature  of  the  reflection  and  using  the  inherent 
properties  to  advantage.  Applications  in  beam  conduction, 
intensity  enhancement,  and  macroscopic  trapping  are  developed 
and  discussed. 

Assumptions 

Several  assumptions  are  used  in  this  thesis.  First, 
only  slow  (~lOOOcm/s)  beams  of  atoms  are  considered.  This 
is  explained  under  the  Basic  Mechanism  section  and  is  due  to 
the  weak  nature  of  the  reflection.  Beam  chopping  is  envi¬ 
sioned  to  eliminate  the  unwanted  high  velocity  portion  of 
the  beam.  Second,  quantum  fluctuations  in  the  force  are 
only  treated  approximately.  A  full  treatment  of  these 
fluctuations  is  beyond  the  scope  of  this  thesis.  Finally, 
atom-atom  interactions  are  ignored  due  to  the  tenuous  nature 
of  the  beams  used  in  current  research.  Further  limited 
assumptions  are  identified  as  they  are  used. 

Scope 

This  thesis  studies  the  ideal  reflection  of  an  atomic 
beam  by  a  near-resonant  evanescent  wave.  Well  collimated 
(.01  radian)  beams  are  considered  throughout.  A  thermal 
velocity  distribution  is  considered  (modified  by  current 
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deceleration  techniques  where  applicable),  or  chopped  to 

discard  the  high  velocity  atoms.  Gravity  is  considered 

throughout  due  to  its  large  effect  on  slow  trajectories. 

Only  physically  realizable  wave  fronts  and  intensities 
o 

('•lwatt/cm  )  are  considered.  Doppler  detuning  of  the  force 
is  considered,  as  well  as  atomic  transitions  out  of  reso¬ 
nance  with  the  laser  beam.  These  factors  limit  the  intensity 
and  length  of  application  of  the  force. 

Summary  of  Current  Knowledge 

At  present  there  are  no  published  works  which  deal 
with  this  exact  problem,  although  there  are  many  in  the 
basic  area  of  resonance-radiation  pressure;  and  "the  subject 
of  atomic  motion  in  resonant  radiation  is  now  a  rapidly 
developing  field  of  research"  (Ref  14).  Dr.  Richard  Cook 
of  the  AFIT  faculty  developed  the  idea  of  using  an  evanes¬ 
cent  wave  in  conjunction  with  the  resonant  pressure  to  create 
an  effective  mirror  for  atoms.  A  good  deal  of  time  was 
spent  rederiving  and  verifying  his  initial  results.  These 
results  are  presented  in  Appendix  A  in  the  form  of  a  joint 
article  which  has  been  accepted  for  publication  in  Optics 
Communications . 

The  theory  of  resonant  radiation  pressure  dates  back 
to  1917  with  a  paper  of  Einstein  titled  "The  Quantum  Theory 
of  Radiation"  (Ref  15).  The  force  was  first  observed  by 
Frish  (Ref  16)  in  1933.  The  advent  of  tunable  lasers 


renewed  interest  in  the  force,  and  much  theoretical  develop¬ 
ment  (Refs  17-24)  and  experimental  verification  (Refs  1-7) 
followed. 

The  latest  interest  in  the  field  has  centered  on 
isotope  separation  (Ref  2),  beam  deceleration  (Ref  7),  and 
the  trapping  of  single  atoms  (Refs  10-13).  This  thesis 
presents  (1)  new  avenues  in  the  last  two  areas 
and  (2)  hope  in  general  beam  control  techniques. 

Approach 

This  thesis  is  an  extension  of  the  results  presented 
in  Appendix  A.  Several  promising  geometries  were  investi¬ 
gated  for  atomic  beams,  reflecting  surfaces  and  incident 
radiation.  The  geometries  which  have  practical  applications 
to  beam  control  or  confinement  are  presented.  The  results 

presented  are  by  no  means  exhaustive  of  the  possible  geome- 

♦  • 

tries,  but  every  effort  was  given  to  finding  the  most 
practical  applications.  Any  lack  of  completeness  is  due  to 
oversight,  as  a  logical  approach  was  taken  in  designating 
promising  geometries. 

The  basic  mechanism  of  evanescent  reflection  is 
examined  first.  This  is  a  straightforward  extension  of 
Appendix  A,  with  an  eye  to  limitations  on  the  force  due  to 
beam  parameters  and  reflector  geometry. 
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An  analysis  of  a  single  reflection  from  a  plane 
surface  is  presented  second.  The  emphasis  here  is  to  develop 
the  single  bounce  velocity  effects  on  atoms  in  the  beam. 

Third,  a  beam  conducting  pipe  is  imagined  and 
examined  as  a  transverse  folding  of  the  plane  surface  around 
the  beam  to  contain  and  channel  the  beam. 

Completely  enfolding  the  beam  and  making  the  reflector 
parabolic  results  in  a  nozzle  which  can  greatly  enhance  the 
beam  intensity  with  little  adverse  effect.  This  is  explained 
in  the  fourth  section  of  the  development. 

Fifth,  the  extension  of  the  transverse  confinement 
of  the  beam  conducting  pipe  to  a  longitudinal  confinement 
leads  to  a  macroscopic  trap  for  slow  atoms  (-^(JOcm/s) . 

This  could  have  direct  application  in  the  field  of  Doppler- 
free  atomic  spectroscopy. 

Finally,  the  results  obtained  are  compared  and  con¬ 
trasted  to  similar  results  for  beam  control  and  trapping  in 
the  current  literature. 

The  BASIC  program  used  to  test  the  thesis  results 
is  presented  in  Appendix  B  to  aid  subsequent  research.  This 
program  follows  the  trajectory  of  an  atom  through  multiple 
bounces  off  three-dimensional  surfaces  which  are  illuminated 
by  very  general  laser  fields. 
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II.  Deve lopment  of  the  Dynamics  of  an 
Electromagnetic  Mirror 

Basic  Mechanism 

The  quantum  mechanical  derivation  of  the  forces 
involved  in  evanescent  reflection  is  given  in  Appendix  A  and 
is  summarized  below.  The  situation  is  as  shown  in  Figure  1. 
However,  note  that  the  axes  are  labeled  differently  in  the 
appendix.  This  is  to  standardize  the  notation  when  gravity 
is  included  in  this  discussion.  The  interaction  is  between 

i 

a  single  freely  moving  atom  and  a  near-resonant  evanescent 
wave. 


Figure  1.  Geometry  of  Basic  Interaction 
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The  atom  is  assumed  to  have  only  two  internal 


quantum  states  and  is  subject  to  no  other  external  forces 

during  the  interaction.  These  two  assumptions  are  valid  for 

most  systems  of  interest.  The  requirement  of  only  two 

quantum  states  can  be  approximated  for  many  atoms  and 

molecules  (Ref  25).  In  a  number  of  cases,  this  requirement 

can  be  achieved  by  the  proper  choice  of  the  transition. 

Another  method  is  to  apply  a  magnetic  field  to  Zeemann  shift 

the  transitions  and  discriminate  against  optical  pumping  of 

the  atom  out  of  resonance.  Results  on  sodium  atoms  (Ref  7) 

have  shown  that  unwanted  transition  rates  can  be  reduced  to 
-4 

~10  times  the  rate  of  the  desired  transition.  This  figure 
will  be  used  to  assess  the  limitations  of  the  force. 

The  neglect  of  noneleetromagnetic  forces  during  the 
interaction  is  very  easy  to  justify.  In  a  normal  interac¬ 
tion,  the  atom  is  affected  by  the  field  for  a  very  short 
time  due  to  the  nature  of  the  evanescent  wave. 

2 

For  an  incident  laser  intensity  of  one  watt /cm  , 
and  a  representative  reflection,  the  atom  undergoes  ~320 
spontaneous  emissions.  This  is  accomplished  in  ~5us.  During 
this  time,  gravitational  acceleration  is  negligible  and  the 
effect  of  a  locally  applied  magnetic  field  for  Zeemann 
shifting  can  be  easily  ignored. 

The  evanescent  wave  is  also  fairly  easy  to  achieve. 

It  arises  from  the  total  internal  reflection  of  a  tunable 
laser  beam  from  a  dielectric  interface.  This  beam  can  be 
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introduced  into  the  dielectric  through  normal  incidence  on 
a  surface  and  directed  toward  the  reflecting  surface.  In 
addition,  the  laser  beam  can  make  numerous  internal  reflec¬ 
tions  by  bouncing  off  dielectric  surfaces.  Consequently,  a 
single  beam  can  be  made  to  illuminate  a  large  area  of  the 
reflecting  surface  with  undiminished  intensity.  This 
arrangement  is  used  in  all  of  the  applications  which  follow. 
The  only  losses  in  this  type  of  system  are  due  to  absorption 
in  the  dielectric,  imperfections  on  the  surfaces,  and  the 
momentum  transfer  to  the  reflected  atoms. 

For  a  plane  electromagnetic  wave  of  frequency  oj 
propagating  in  a  dielectric  of  refractive  index  n,  the 
evanescent  wave  will  have  the  form 


fi(x,t)  =  S£e“azcos{wt-kx} 


where  £  =  the  unit  polarization  vector, 

£  =  wave  amplitude  at  z=0, 
a  =  oa(n2sin20-l)^/c , 
k  =  amsin0/c, 

and  0  =  angle  of  incidence. 

The  interaction  consists  of  two  forces.  The  first 

i 

gives  rise  to  the  zero-order  effect:  the  reflection.  This 
is  the  force  normal  to  the  surface,  and  is  given  by 


2aR(A-kv„)jr(z) 


4(A-kv .)  +A  +2fl  (z) 


where  A  =  0)-a>o  (detuning) 

ft  =  y£exp(-az)/fi  (Rabi  frequency) 

2  3  3 

A  =  4y  o)q/3Rc  (Einstein  A  coefficient) 

u0=  atomic  transition  frequency 
and  y  =  dipole  transition  moment. 

Note  that  for  reflection  this  force  must  be  positive,  neces¬ 
sitating  a  positive  A.  That  is,  the  laser  must  be  tuned 
above  the  transition  frequency  of  the  atom. 

The  second  force  is  parallel  to  the  surface  and  in 
the  direction  of  the  evanescent  wave  vector.  This  force  is 
given  by 

F  =  RkAft2(z) 

x  4(A-kv  )2+A2+2ft2(z)  (3) 

X 

This  is  the  first-order  effect  of  the  reflection,  and  serves 
to  give  an  added  "kick"  to  the  atom  in  the  direction  of  the 
evanescent  wave  vector.  It  is  this  force  which  can  be  used 
to  decelerate  or  redirect  parts  of  an  atomic  beam.  It  is  a 
hinderance,  however,  in  collimation  and  trapping. 

Two  important  results  follow  from  Eqs  (2)  and  (3). 
The  first  gives  the  maximum  normal  velocity  which  can  be 
reflected.  This  is  the  major  limitation  on  the  force  and  is 
given  by 


max  _  rRA'  nM.  2y"£~ 

—  t  u  ln(l+  o  22 

z  M  a  [4(A-kv  )2+A  J 
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where  M  is  the  mass  of  the  atom,  and  L'  includes  doppler 
effects.  The  second  result  directly  relates  the  strength  of 
the  kick  to  the  strength  of  the  reflection.  The  ratio  of 
these  two  forces  is  designated  R,  and  is  given  by 


p  _  _x  _  kA 
K  Fz  2cxA 


nAsinO 

2A>(n2sin29-l)i 


This  ratio  is  independent  of  position  which  means  that  the 

atom  will  always  experience  a  kick  which  is  proportional  to 

its  initial  velocity  normal  to  the  surface. 

To  get  a  feel  for  these  two  results,  consider  the 

2  2 

case  of  a  sodium  atom  utilizing  the  standard  3  P3/2 

transition.  Then  M  =  3.817xl0~23g  ,  <i)Q  =  5 . 090xl014s-1  , 

7  —1  — lfi 

A  =  6.22x10  s'  ,  and  p  =  1.003x10  .  Using  a  nominal 

o 

laser  irradiance  I  *  lwatt/cm  ,  the  field  strength  £,  and 
the  Rabi  flopping  freuqency  0  depend  only  on  the  refractive 
index  n  and  the  angle  of  incidence  0  of  the  laser.  Plots  of 


C2  =  128Trn2cos20I/(n2-l)(n+l)2c 


no  =  vz/h 


are  shown  in  Figure  2  for  various  values  of  n  and  0  (0Q  is 
the  Rabi  frequency  at  z=0).  For  the  purposes  of  this  thesis, 
values  of  n=1.5  and  0=45°  give  a  good  trade-off  between 
available  dielectric  media  and  low  angles  of  incidence.  A 
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i 


low  angle  of  incidence  is  desirable  for  multiple  reflections 
as  outlined  earlier. 

With  n  and  9  fixed,  the  maximum  normal  velocity  and 
the  kick  ratio  depend  only  on  the  detuning  A  (plus  doppler 
effects).  Plots  of  this  dependence  are  given  in  Figure  3. 
The  figure  is  split  into  low  and  high  detuning  regimes  to 
show  the  ranges  for  maximum  kick  and  reflection.  Because 
Doppler  effects  must  be  considered  in  the  detuning,  uncer¬ 
tainties  are  shown  for  both  regimes  appropriate  to  the 

expected  range  of  longitudinal  velocities  (after  velocity 

8  —1 

selection).  Nominal  detuning  values  of  5x10  s  for  a 

9  -1 

moderate  kick  and  5x10  s  for  reflection  of  maximum  normal 
velocity  will  be  used  henceforth. 

Plane  Mirror  Analysis 

Semi-Classical  Effect.  The  easiest  case  of  electro¬ 
magnetic  reflection  to  analyze  is  that  of  simple  reflection 
from  a  plane  surface.  Therefore,  consider  a  dielectric 
interface  illuminated  from  within  by  a  plane  wave  as  shown 
in  Figure  1.  Atoms  in  the  beam  which  have  a  normal  velocity 
component  less  than  the  maximum  (uraax)  dictated  by  I,  n,  9, 
and  A  will  be  reflected  with  an  equal  (but  opposite)  normal 
component.  During  this  reflection,  however,  they  will  gain 
(or  lose)  momentum  from  the  spontaneous  interaction  with  the 
field.  This  momentum  will  be  gained  parallel  to  the  inter¬ 
face  and  in  the  direction  of  the  k-vector  of  the  evanescent 
wave. 
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v  max  and  R  for  Sodium 
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For  counter-directed  radiation,  this  will  slow  all 
atoms  which  are  reflected.  For  an  ideal  beam  of  thermal 
velocity  distribution  and  grazing  incidence  at  angle  <f>,  the 
results  are  as  follows.  All  atoms  with  velocity  less  than 
umax/sin<j>  will  be  reflected.  Atoms  with  velocities  greater 
than  this  value  will  strike  the  dielectric  interface  and 
either  be  absorbed,  scattered  elastically,  or  inelastically . 

Atoms  which  are  absorbed  or  inelastically  reflected 
will  damage  the  surface.  Elastically  scattered  atoms  will 
lose  a  momentum  increment  proportional  to  the  time  they 
interact  with  the  field.  This  is  a  complex  function  roughly 
inversely  proportional . to  the  atom's  initial  velocity.  No 
reliable  estimate  of  the  number  of  elastically  scattered 
atoms  is  possible;  therefore,  atoms  with  normal  velocities 
greater  than  umax  will  be  omitted  from  further  consideration, 
except  as  a  limit  to  thd  interaction. 

Atoms  with  velocities  less  than  umax  will  be  reflected 
by  the  evanescent  wave.  In  addition  to  a  reversal  of  their 
initial  normal  velocity,  these  atoms  will  lose  longitudinal 
velocity  due  to  spontaneous  emissions.  The  velocity  incre¬ 
ment  lost  is  equal  to  2R  times  their  initial  normal  velocity, 
i.e.  , 


<$v  =  -2Rvsin<f>  (8) 

This  effect  can  be  used  to  enhance  the  low  velocity  portion 
of  a  thermal  beam,  by  using  grazing  incidence  from  an 
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interface.  The  population  increase  due  to  the  downward 
shifting  of  velocities  from  a  thermal  beam  is  shown  in 
Figure  4.  The  result  is  a  22%  increase  in  the  lew  velocity 
population  under  nominal  conditions. 

In  general,  for  laser  illumination  at  an  angle  y 
from  the  atomic  beam  direction  (measured  in  the  plane  of  the 
interface),  the  velocity  increment  will  be 

<Svx  =  2Rvcosysin<}>  (9) 

The  cases  of  most  interest,  however,  are  counter-propagation 
(y  =  180°)  for  beam  deceleration,  co-propagation  (y  =  0)  for 
acceleration,  and  cross  propagation  (y  =  90°)  for  transverse 
deflection.  Other  angles  will  yield  a  combination  of  accel¬ 
eration  and  deflection.  In  the  beam-conducting  pipe  and 
trapping  applications,  the  effects  are  resolved  into  longi¬ 
tudinal  and  transverse  components  to  produce  elastic  bounces 
for  confinement  and  longitudinal  kicks  to  accelerate  the 
beam.  For  the  enhanced  nozzle  application,  the  angle  is 
unimportant  except  for  practical  considerations  of  laser 
geometry.  This  is  due  to  the  suppression  of  the  kick  for 
this  application. 

Quantum  Considerations.  One  effect  which  has  been 
overlooked  thus  far  is  a  velocity  spreading  due  to  the  nature 
of  the  force.  The  kick  arises  from  spontaneous  decays  of 
the  atoms  which  occur  in  random  directions.  The  atom 
experiences  a  net  kick  in  the  direction  of  the  laser  because 
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it  acquires  the  momentum  of  one  photon  on  each  absorption. 
It  then  loses  this  momentum  in  a  random  direction  resulting 
in  an  average  momentum  gain  in  the  direction  of  the  laser's 
k-vector.  The  random  nature  of  the  emission  produces  a 
statistical  gaussian  spreading  of  the  atom's  velocity. 

An  evaluation  of  the  results  of  this  spreading 
requires  a  computation  of  the  numbers  involved.  On  each 
photon  absorption,  the  atom  gains  an  increment  of  velocity 
given  by 


6v  =  liw/Mc  (10) 

The  number  of  spontaneous  emissions  N  and  the  time  required 
for  these  emissions  (essentially  the  interaction  time)  fit 
are  then  given  by 


N  =  2RvzMc/Bu>  (11) 

and 

fit  =  N/A  =  2RvzMc/fiu)A  (12) 

The  distribution  in  velocities  is  then  characterized  by  a, 
the  room-mean-square  deviation,  given  by  (Ref  26) 

a  =  /fifiv  (13) 

For  the  two  cases  of  interest,  moderate  kick  and 
maximum  reflection,  the  velocity  increment  is  the  same: 
fiv  =  0.47cm/s  per  emission.  For  the  moderate  kick  case  and 
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the  limiting  normal  velocity,  N  =  320  emissions, 

St  =  5.2ys  ,  and  a  =  8.4cm/s  .  For  the  case  of  maximum 
reflection  (again  at  limiting  normal  velocity),  N  =  64 
emissions,  St  =  lys  ,  and  a  =  3.75cm/s  . 

This  derivation  is  only  approximately  correct.  A 
more  complete  treatment  (Ref  14)  reveals  that  the  induced 
spreading  is  larger  than  expected.  For  the  mechanism  of 
this  thesis,  there  is  an  increased  spreading  in  the  direc¬ 
tion  of  propagation  of  the  evanescent  wave.  There  is  also 
an  induced  spread  (or  diffusion)  normal  to  the  surface  due 
to  the  gradient  in  the  Rabi  frequency.  A  careful  analysis 
of  these  spreadings  was  done  using  the  diffusion  coefficients 
obtained  by  Cook  (Ref  14).  For  the  geometry  of  interest, 
this  required  numerical  integration  of  the  various  diffusion 
coefficients  with  respect  to  time  for  a  path  in  the  evanescent 
potential  well. 

The  spontaneous  root-mean-square  deviations  vary  as 
the  square  root  of  the  initial  normal  velocity  of  the  atom. 

The  induced  deviation  in  the  direction  of  propagation  is 
also  proportional  to  the  square  root  of  the  initial  normal 
velocity,  but  is  ~25%  larger  than  the  spontaneous  deviation. 
The  induced  deviation  normal  to  the  surface  is  composed  of 
a  linear  and  a  cubic  term  in  the  initial  velocity.  There¬ 
fore  the  induced  diffusion  normal  to  the  surface  is  much 
more  important  at  higher  atomic  velocities!  The  coefficients 
for  spontaneous  and  induced  diffusion  are  given  in  Appendix  B 
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in  the  BOUNCE  routine  listing  for  the  two  nominal  detunings 
of  this  thesis.  For  other  detunings,  the  coefficients  should 
scale  roughly  as  the  inverse  of  the  fourth  root  of  the 
detuning . 

Because  the  induced  diffusion  normal  to  the  surface 
varies  as  the  cube  of  the  initial  velocity,  this  term  sets 
an  effective  upper  bound  on  the  initial  velocities.  For  the 
moderate  kick  detuning,  this  induced  diffusion  is  larger 
than  the  spontaneous  diffusion  for  velocities  above  ~70cm/s. 
For  this  detuning,  the  root-mean-square  deviation  is  ~12cm/s 
for  an  initial  velocity  of  lOOcm/s.  For  the  maximum  reflec¬ 
tion  detuning,  the  cross-over  is  at  90cm/s  and  the  deviation 
at  lOOcm/s  is  only  1.8cm/s.  For  this  larger  detuning,  the 
deviation  does  not  exceed  lOcm/s  for  initial  velocities 
below  175cm/s. 

The  deviation  in  the  propagation  direction  is  modest 
but  will  tend  to  enhance  the  low  velocity  portion  of  the  beam 
due  to  the  preponderance  of  higher  velocities  in  the  Maxwell- 
Boltzraann  distribution.  The  effect  on  the  population  is 
shown  in  Figure  5  for  the  moderate  kick  case.  This  figure 
is  only  an  approximation  as  the  integrals  involved  are 
analytically  unsolvable  and  numerically  unwieldy. 

There  are  two  other  effects  of  this  gaussian  velocity 
spreading  which  should  be  mentioned;  spreading  of  the  trans¬ 
verse  and  normal  velocities.  The  spreading  transverse  "~o 
the  atomic  beam  axis  will  tend  to  disperse  the  beam.  This 
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dispersion  will  be  small,  on  the  order  of  inherent  disper¬ 
sions,  even  for  very  slow  beams.  In  confinement  schemes, 
it  will  tend  to  heat  the  atoms;  but  the  effect  should  be 
minimal.  Spreading  of  the  velocity  normal  to  the  interface 
will  have  two  effects.  The  first  effect  is  the  spread  of 
reflected  velocities  just  as  in  the  other  directions.  The 
second  effect  arises  from  velocity  fluctuations  during  the 
reflection.  Because  of  these  fluctuations,  the  atom  loses 
the  "memory"  of  its  initial  normal  velocity.  This  will 
complicate  the  reflection  process  and  change  the  kick  and 
the  diffusion  in  all  directions.  Fortunately,  this  effect 
is  very  small  and  will  be  ignored. 

Beam  Conducting  Pipe 

Longitudinal  Characteristics.  A  straightforward 
extension  of  the  plane  mirror  in  the  last  section  will  pro¬ 
duce  a  surface  which  can  be  used  to  conduct  a  slow  atomic 
beam  from  one  place  to  another.  For  this  purpose,  gravita¬ 
tional  trajectories  must  be  considered  between  interactions. 
The  present  section  considers  only  ideal  beams  of  atoms, 
perfectly  collimated,  and  interacting  with  the  evanescent 
wave  only  during  an  incremental  bounce.  The  first  condition 
will  be  relaxed  to  encompass  dispersive  beams  in  a  later 
section.  The  second  condition  is  essentially  true  for 
realizable  evanescent  waves  with  n  =  1.5  and  0  =  45°  . 
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3  -1 

These  give  a  =  6x10  cm  ,  and  the  interaction  region  is 
on  the  order  of  2  microns  deep. 

Consider  then  a  horizontal  mirror,  illuminated  as 
described,  and  an  ideal  atomic  beam  issuing  horizontally 
from  a  nozzle  a  distance  h  above  the  mirror.  The  atoms 
will  then  fall  due  to  gravity  (g)  and  strike  the  mirror 
surface  at  a  distance  1  given  by 

1  =  v/2h/g  (14) 

The  atoms  will  be  reflected  at  this  point  for  h  such  that 
the  normal  velocity  (due  to  gravity)  is  less  than  umax  as 
defined  earlier.  This  limits  h  to 

h  <  umax2/2g  (15) 

For  umax  ~  300cm/s  ,  this  gives  a  maximum  h  of  ~45cm.  This, 
however,  results  in  excessive  values  of  1.  For  a  maximum 
v  =  lOOOcm/s  impacting  at  1  =  lm  ,  h  becomes  5cm.  This 
gives  a  lower  normal  velocity  and  thereby  a  lesser  decelera¬ 
tion  on  each  bounce.  A  tradeoff  must  be  made  based  on  the 
initial  velocity  composition  of  the  beam,  the  available 
mirror  length,  and  the  purpose  of  the  deceleration.  Low 
velocities  are  of  the  most  interest  for  applications  in 
atomic  clocks,  ultrahigh  resolution  spectroscopy,  etc. 

Each  atom  would  be  decelerated  by  the  same  amount  on 
each  bounce : 


<5v  =  2R/5gE 


(16) 


For  a  mirror  of  length  L,  this  results  in  a  shift  of  the 
initial  velocity  distribution  to  lower  velocities.  Indeed 
for  initial  velocities  below  a  cutoff  determined  by  L,  the 
atoms  will  be  retro-ref lected  toward  the  nozzle.  This  shift 
is  illustrated  in  Figure  6.  Note  that  the  shift  to  lower 
velocities  is  incremental.  In  addition,  the  atoms  are 
smeared  between  z  =  0  and  z  =  h  based  on  the  position  of 
the  last  bounce  before  they  reach  the  end  of  the  mirror. 

This  vertical  distribution,  coupled  with  the  incremental 
nature  of  the  longitudinal  velocities  at  a  given  height 
could  be  used  to  select  a  very  small  velocity  segment  of  the 
beam. 

An  interesting  result  is  obtained  by  solving  Eqs 
(14)  and  (16)  for  the  time  average  of  the  acceleration 
experienced  by  an  atom.  This  yields 

a  =  ( 5v)v/21  =  -Rg  (17) 

independent  of  h  and  v.  That  is,  the  reflections  act  as  a 
simple  (but  incremental)  deceleration  of  the  beam.  By 
tilting  the  surface  (and  nozzle)  downward  this  deceleration 
can  be  counteracted  by  gravity  to  produce  an  acceleration- 
free  conduction  of  the  atomic  beam.  The  same  is  true  for 
co-propagation  of  the  laser  and  an  upward  tilt  of  the  sur¬ 
face.  Because  the  geometry  separates  into  motion  along  the 
surface  and  a  gravitational  displacement,  the  angle  of 
tilt  <f>  for  acceleration-free  conduction  is  simply 
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Figure  6.  Velocity  Shift  from  Beam  Pipe.  For  each  reflec¬ 
tion  a  population  group  is  shifted  by  ~50cm/s.  This  results 
in  gaps  at  some  velocities  and  greatly  increased  population 
at  others. 
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<t>  =  tan_1R 


(18) 


For  n  =  1.5  ,  0  =  45°  ,  and  A  =  5xl08s-^  ,  this  yields  a 
tilt  angle  <j>  =  10.6°  . 

By  not  tilting  the  surface,  the  incremental  decele¬ 
ration  can  be  used  to  increase  the  low  velocity  portion  of 
the  beam  by  shifting  the  more  populous  segments  of  the  beam 
to  lower  velocities.  The  advantage  of  this  method  over 
deceleration  by  gravity  is  the  elimination  of  unwanted 
velocities  at  the  point  of  observation  (by  the  incremental 
vertical  shifting),  coupled  with  the  transverse  confinement 
of  the  beam  to  be  discussed  in  the  next  section. 

Figure  6  shows  the  dramatic  improvement  in  the  low 
velocity  population  due  to  the  incremental  shift.  Although 
the  pipe  lengths  portrayed  were  picked  for  convenience 
rather  than  population  enhancement,  they  give  rise  to  signif 
icant  enhancement.  For  a  pipe  length  of  two  meters,  the 
total  population  of  atoms  with  longitudinal  velocities 
below  lOOcm/s  is  increased  to  over  23  times  its  original 
value.  For  the  four  meter  pipe,  the  improvement  is  35  times 
the  original  population.  By  a  proper  choice  of  pipe  length 
coupled  with  velocity  chopping  at  the  output  end,  a  large 
sample  of  atoms  with  very  well  defined  velocity  spread  can 
be  produced. 

The  gaussian  velocity  spread,  due  to  the  quantum 
nature  of  the  force,  produces  a  beneficial  result.  The 
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spreading  will  tend  to  round  the  sharp  edges  of  the  shifted 
distributions  shown  in  Figure  6.  However,  the  fact  that  the 
initial  population  is  skewed  toward  the  high  velocity  regime 
will  result  in  further  enhancement  of  the  low  velocity  por¬ 
tion.  The  amount  of  this  enhancement  is  dependent  on  the 
pipe  length  used  and  no  estimate  was  made. 

Transverse  Characteristics.  Consider  now  the  motion 
of  atoms  transverse  to  the  laser  direction  in  the  above  dis¬ 
cussion.  Because  the  induced  kick  is  in  the  direction  of 
the  evanescent  wave's  k-vector,  transverse  motion  is  unaf¬ 
fected  except  for  the  gaussian  spreading  due  to  the  nature 
of  the  kick  as  described  above.  The  effect  of  this  trans¬ 
verse  velocity  spreading  will  be  ignored  for  the  moment  and 
discussed  at  the  end  of  this  section. 

The  problem  arises  as  to  how  the  mirror  surface  can 
be  shaped  perpendicular  to  the  atomic  beam  axis  in  order  to 
contain  and  conduct  the  beam.  Because  these  bounces  do  not 
change  the  transverse  speed  of  the  atoms,  it  becomes  a 
problem  in  standard  kinetics. 

For  an  atom  with  transverse  velocity  vy  starting  at 
a  height  h,  its  velocity  components  at  any  later  position 
(y,z)  are  given  by 

=  /2g(h-z ) 


v 


(19) 


Vy  =  -yg//2g(h-z) 


(20) 


For  retroref lection  in  the  transverse  direction,  the  surface 
normal  must  match  the  velocity  vector.  This  determines  z  as 
a  function  of  y  for  the  surface.  The  result  is 

z  =  h  -  /h2  -  y2/2  (21) 


a  simple  ellipse. 

This  surface  will  return  all  atoms  to  their  original 
apex  at  y  =  0  and  z  =  h  for  transverse  velocities  up  to 
the  maximum  limit  for  reflection.  For  the  moderate  kick 
arrangement  (umax  =  289cm/s)  and  normal  collimation  (.01 
radian)  this  means  that  beams  chopped  for  velocities  above 
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2.9x10  can  be  conducted  without  loss  due  to  lateral  disper¬ 
sion.  This  analysis,  however,  ignores  the  effect  of  the 
vertical  component  of  the  dispersion.  Simulations  using  the 
BASIC  program  in  Appendix  B  show  that  the  transverse 
channeling  begins  to  weaken  for  vertical  components  at  the 
beam  axis  in  excess  of  lOcm/s.  This  weakening  is  manifested 
as  a  shift  in  the  apex  position  away  from  y  =  0  .  This 
shift  results  in  a  mismatch  between  the  atom's  velocity 
components  and  the  normal  vector  of  the  surface.  The  result 
is  a  coupling  between  the  y  and  z  oscillations  of  the  tra¬ 
jectory.  The  position  of  the  apex  shifts  after  each  bounce 
in  an  unstable  pattern.  The  observed  tendency  is  toward 


higher  z  values  initially,  followed  by  a  decay  toward  lower 
values.  This  pattern  repeats  itself  on  the  order  of  four 
bounces . 

In  order  to  examine  and  control  this  oscillation  an 
attempt  was  made  to  design  a  surface  which  would  compress 
the  lateral  motion  of  the  beam.  This  would  effectively  focus 
the  beam  laterally  and  greatly  enhance  the  intensity  avail¬ 
able  for  subsequent  observation.  The  approach  used  was  a 
modification  of  Eq  (21)  to  incorporate  a  squeezing  parameter 
8.  This  results  in 

z  =  h  -  /h2  -  8y2/2  (22) 

Results  from  this  approach  were  disappointing.  For  8  <  1  , 
the  oscillations  were  intensified.  For  8  >  1  ,  the  oscil¬ 
lations  initially  reversed,  tending  toward  lower  apexes  but 
then  continued  unabated.  These  results  are  included  to  pre¬ 
vent  repetition  of  these  fruitless  efforts. 

Turning  to  the  gaussian  velocity  spreading  of  the 
beam  due  to  the  spontaneous  nature  of  the  force,  the  effect 
on  this  transverse  confinement  scheme  is  minimal.  On  each 
bounce  all  velocity  components  will  acquire  a  deviation 
based  on  the  reflection  parameters  used. 

The  result  of  the  spreading  of  the  z  component  is  a 
spreading  of  the  apex  height,  and  consequent  mismatch  with 
the  design  parameters  of  the  pipe.  For  a  vertical  velocity 
v  and  change  in  vertical  velocity  Sv,  the  resultant  height 


change  is  given  by 


5h  =  (2v  6v  +  5 v2 ) /2g  (23) 

z 

The  dependence  of  v  on  v  is  more  complicated,  but  is 

z  y 

characterized  by 


2  4 

v  v 

v  =  ~/2gh  +  — y- — - 3  +  •  •  • 

z  4/2gh  32(  /2gh;  (24) 

which  leads  to  a  height  change  of 

6h  *  6v2/2g  +  . 84 Jh8v/g  -  .84v26v//fr  (25) 

For  the  beam  pipe  parameters  described  and  moderate 
kick  detuning,  lateral  velocities  up  to  275cm/s  will  be 
reflected.  The  spreading  of  the  vertical  velocities  will 
result  in  a  spread  in  the  apex  height  on  the  order  of  0.6cm. 
This  is  well  within  the  normal  spread  of  heights  due  to  the 
dispersion  of  the  beam  for  beam  speeds  ~700cm/s  contemplated 
in  this  application. 

The  spread  of  velocities  in  the  y  direction  is  some¬ 
what  more  serious.  It  results  in  a  transverse  shift  in  the 
apex  position  and  consequent  mismatch  with  the  pipe.  Com¬ 
puter  simulation  has  shown,  however,  that  the  pipe  as 
described  can  tolerate  transverse  velocity  shifts  on  the 
order  of  lOcm/s  which  is  roughly  the  anticipated  sigma  for 
lOOcm/s  transverse  velocities.  This  means  that  68%  of  the 


deviations  can  be  handled  by  the  transverse  shape  of  the 
pipe. 


Next,  consider  the  effect  of  beam  dispersion  on  the 

use  of  such  a  beam  pipe.  The  design  is  very  sensitive  to 

vertical  velocity  on-axis.  This  limitation  is  '•lOcm/s, 

which  for  collimation  of  0.01  radian  limits  the  top  beam 
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speed  to  10  cm/s.  The  lateral  requirements  are  not  as 
severe.  The  pipe  will  contain  lateral  velocities  up  to 
lOOcm/s  with  little  adverse  effect.  This  can  be  exploited 
by  using  a  wide  thin  opening  at  the  oven  coupled  with  a  small 
opening  at  the  pipe.  This  gives  a  fan  shaped  beam  which  is 
vertically  narrow.  This  will  result  in  a  lOfold  increase  in 
beam  intensity  over  standard  collimation,  although  the  out¬ 
put  would  still  contain  the  lateral  velocity  spread.  How¬ 
ever,  this  will  allow  placement  of  detectors  in  a  very  pure 
longitudinal  velocity  segment  of  the  beam. 

One  problem  which  has  not  been  addressed  thus  far  is 
the  matching  of  the  laser  radiation  to  the  geometry  of  the 
beam  conducting  pipe.  Ignoring  the  lateral  shaping,  align¬ 
ment  is  simple.  The  mirror  is  a  flat  plane  which  becomes 
the  top  surface  of  a  thin  (~lcm)  dielectric  sheet.  The 
beam  output  end  of  the  sheet  is  then  given  a  45°  bevel  and 
the  laser  is  introduced  to  the  sheet  normal  to  this  bevel. 
When  the  sheet  is  curved  around  the  beam  ax^  to  provide 
containment,  the  situation  is  not  quite  so  simple. 


The  inner  surface  of  the  pipe  is  an  ellipse,  as 
described  earlier.  The  problem  is  to  find  the  proper  outer 
shape  for  the  pipe  to  match  the  inner  curvature.  The  out¬ 
put  end  could  then  be  beveled  as  before,  and  the  laser 
expanded  and  focused  to  strike  the  normal  of  the  bevel.  The 
radiation  would  then  travel  down  the  pipe  much  as  before, 
focusing  in  toward  the  inner  surface,  reflecting,  expanding 
toward  the  outer  surface  for  another  reflection,  and 
repeating  this  process  as  it  traverses  the  length  of  the 
pipe. 

The  proper  shape  for  the  outer  surface  is  another 
ellipse  centered  on  the  beam  axis.  This  ellipse  is  a 
distance  d  larger  than  the  inner  surface  at  all  points.  The 
resulting  equation  is 

(z-h)2  «  (h+d)2-x2(h+d)2/(/Zh+d)2  (26) 

The  laser  irradiance  varies  as  1/d,  but  disturbing 
waveguide  effects  occur  for  d  less  than  ~lmm.  Therefore,  a 
determination  of  the  size  of  d  must  include  consideration 
of  the  laser  resources  available  and  the  user's  end  purpose 
for  the  pipe. 

Enhanced  Nozzle 

One  application  of  resonant  reflection  which  could 
be  of  great  help  in  the  study  of  slow  atoms  is  a  reflective 
nozzle  to  aid  beam  collimation.  Because  of  the  extremely 


low  populations  of  slow  atoms  in  a  thermal  beam,  any 
improvement  in  collimation  without  concomitant  losses  would 
greatly  increase  the  number  of  atoms  for  study. 

By  forming  a  parabolic  reflector  around  the  beam 
with  a  suitable  dielectric,  such  a  nozzle  can  be  constructed 
The  arrangement  is  shown  in  Figure  7. 


Figure  7.  Nozzle  Geometry.  A  beam  of  atoms 
diverges  from  the  slit  and  is  reflected  from  a  dielectric 
paraboloid  illuminated  from  within. 

For  this  application,  the  induced  kick  on  each 
bounce  is  a  hindrance.  Therefore,  use  of  a  large  detuning 
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(5x10  s  )  for  maximum  reflection  and  minimum  kick  is  antic¬ 
ipated.  This  limits  the  maximum  normal  velocity  at  the 
nozzle  surface  to  455cm/s.  For  a  small  beam  source  at  the 
focus  of  the  parabola  (x=p)  [see  Figure  7] ,  and  atoms 
issuing  at  angle  9;  the  distance  to  reflection  (r),  and  the 
x  and  y  coordinates  of  reflection  are  given  by 


2p(cos(9 )+l/sin2( 9) 

(27) 

x  =  p  +  rcos(0) 

(28) 

y  =  rsin  (0) 

(29) 

At  angle  0,  the  maximum  longitudinal  velocity  which  can  be 
reflected  is  given  by 

v™ax  =  455/tan( 0 ) sin (9/2)  (30) 

Thus  x  velocities  up  to  3000cm/s  can  be  reflected  for  angles 
up  to  30°. 

In  any  application  of  this  resonant  nozzle  a  trade¬ 
off  must  be  made  between  the  maximum  longitudinal  velocity 
and  the  maximum  angle,  and  between  the  length  of  the  nozzle 
and  the  minimum  angle.  To  explore  the  utility  of  the  nozzle, 
nominal  values  of  8  and  30  degrees  were  chosen  for  the 
minimum  and  maximum  angles  with  p  =  0.1cm  . 

These  parameters  give  a  maximum  longitudinal  velocity 
of  3000cm/s  as  described  above,  a  nozzle  length  of  20.5cm 
and  an  output  radius  of  2.9cm.  This  radius  is  an  order  of 
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magnitude  less  than  that  for  a  collimated  beam  of  3000cm/s 
atoms  at  lm.  In  addition,  the  major  portion  of  the  atoms 
are  reflected  from  the  larger  angles  (where  the  radial  dis¬ 
placement  is  the  least).  For  angles  above  20  degrees,  y  is 
less  than  1.1cm.  At  the  output  end  of  the  nozzle,  ignoring 
the  induced  kick  which  will  be  discussed  later,  the  atoms 
are  perfectly  collimated  and  ready  for  observation  or  con¬ 
duction  via  a  beam  pipe  for  deceleration,  etc.  The  most 
startling  change  in  the  beam  is  the  enhancement  of  the  popu¬ 
lation  available  for  further  manipulation.  By  collecting 
and  collimating  the  atoms  at  angles  between  8  and  30  degrees, 
the  number  available  is  2484  times  that  contained  in  the 
normal  .01  radian  sample!  By  using  an  elliptic  (rather  than 
parabolic)  nozzle,  the  input  beam  can  be  focused  onto  a 
small  spot  for  filling  an  atomic  trap  or  further  manipula¬ 
tion  and  detection.  * 

Note  that  gravitational  displacement  has  not  been 
mentioned  in  the  above  discussion.  This  is  due  to  the 
anticipated  velocity  range  used  in  this  application.  For  a 
longitudinal  velocity  of  700cm/s,  the  atom  will  only  drop 
0.42cm  in  the  20.5cm  length  of  the  nozzle.  The  drop  for 
higher  velocities  trails  off  as  the  inverse  square  of  the 
velocity.  This  effect  will  be  negligible  when  compared  with 
the  focusing  effects  of  the  nozzle  and  the  quantum  fluctua¬ 
tions  due  to  the  nature  of  the  force.  The  requirement  for 
moderate  initial  velocities  can  be  negated  by  using  a  beam 


pipe  for  deceleration  at  the  output  end  of  the  nozzle. 

Changes  in  the  nozzle  parameters  will  also  reduce  the  gravi¬ 
tational  effect.  By  increasing  the  minimum  angle  to  10 
degrees,  the  nozzle  length  is  decreased  to  13.1cm  from  20.5 
with  reduction  of  the  population  multiplier  to  2375  from 
2484.  It  appears  that  a  more  serious  limitation  is  the  proper 
construction  and  alignment  of  the  nozzle  and  the  beam 
origin . 

Returning  to  the  kick  induced  on  the  bounce,  there 
are  three  effects  worthy  of  note.  Assuming  counter  propa¬ 
gation  of  the  laser,  these  are  a  deceleration  of  the  atoms, 
an  over-focusing  due  to  the  angle  of  reflection,  and  the 
inevitable  qaussian  velocity  spreading  in  all  directions. 

The  deceleration  is  given  by 

6v  =  2RvYtan8sin(8/2)cos(8/2)  (31) 

For  the  parameters  above,  this  is  a  fairly  modest  effect. 

For  vx  =  3000cra/s  and  8=30°  ,  the  deceleration  is  only 
16cm/s,  decreasing  rapidly  for  lower  velocities  and  angles. 

The  kick  also  has  a  component  in  the  radial  direc¬ 
tion.  This  serves  to  increase  the  focusing  effect  of  the 
parabola,  and  is  given  by 

<5Vj  =  2Rvxtan8sin^(8/2)  (32) 

For  the  nominal  parameters  above,  this  is  never  more  than 
4.4cm/s,  and  is  typically  less  than  lcm/s.  This  is  well 


within  the  allowable  range  of  deviations  and  less  than  the 
effect  of  the  gaussian  velocity  spread. 

The  effect  of  the  gaussian  velocity  spread  can  be 
resolved  into  three  components.  For  the  maximum  kick 
detuning  used  in  this  application,  the  spread  is  character¬ 
ized  by  a  root-mean-square  deviation  of  lOcm/s  times  the 
cube  of  the  ratio  of  the  normal  velocity  over  175cm/s. 

The  spread  in  the  axial  direction  can  be  ignored  relative 
to  the  axial  velocities  ranging  in  the  thousands.  There  is, 
of  course,  the  beneficial  effect  on  the  lower  velocity 
portion  due  to  the  skewed  initial  population.  The  spreading 
of  velocities  in  the  <p  direction  can  also  be  ignored  as  it 
is  negligible  compared  to  construction  and  misalignment 
effects.  The  radial  component  of  the  spreading  is  alone 
worthy  of  further  examination. 

The  radial  spreading  component  is  the  major  innate 
cause  of  defocuslng,  second  only  to  misalignment  of  the  beam 
source.  The  effect  is  only  significant  for  longitudinal 
velocities  above  ~750cra/s.  For  this  longitudinal  velocity, 
the  radial  velocity  induced  at  two  sigma  (which  encompasses 
95%  of  the  deviations)  is  only  7.5cm/s,  within  allowable 
limits  for  a  subsequent  beam  pipe  and  corresponding  to  a 
0.1  radian  collimation  at  750cm/s. 

The  problem  of  laser  illumination  is  more  serious. 
The  nozzle  should  be  formed  with  a  parabolic  interior,  and 
an  outer  shape  which  ensures  a  constant  angle  between  the 


laser,  multiply  reflecting  down  the  dielectric  interior,  and 
the  interior  interface.  The  equation  for  the  outer  shape 
unfortunately  involves  an  integral  with  no  analytic  solution 
In  fact,  there  is  no  outer  shape  which  will  satisy  the 
requirement  exactly.  The  best  solution  is  a  thin-walled 
paraboloid  which  will  allow  very  little  variation  in  the 
interface  angle  over  the  length  of  the  pipe.  As  long  as  the 
wall  is  thick  enough  to  avoid  spill  over  problems  associated 
with  thin  waveguides,  the  variation  should  not  be  signifi¬ 
cant.  One  factor  working  in  the  designer's  favor  is  the 
inherent  narrowing  of  the  paraboloid.  This  insures  that  the 
effective  illumination  will  be  more  intense  near  the  base 
of  the  paraboloid  due  to  the  smaller  effective  area.  This 
will  help  compensate  for  the  varying  angles  of  incidence. 

For  the  parameters  above  and  a  0.25cm  thick  wall,  no  insur¬ 
mountable  problems  are  foreseen. 


Trap  for  Slow  Atoms 

The  lateral  confinement  scheme  used  in  the  beam  con¬ 
ducting  pipe  is  the  starting  point  for  a  scheme  to  trap  and 
hold  slow  atoms  for  observation.  The  system  is  already 
bound  in  four  directions:  plus  z  by  gravity,  minus  z  by  the 
reflecting  surface,  and  plus  and  minus  y  by  the  elliptical 
walls  of  the  cavity.  The  remaining  problems  are  in  over¬ 
coming  the  induced  kick,  designing  the  x  variation  of  the 
trap,  evaluating  the  effectiveness  and  usefulness  of  the 


trap,  and  finding  methods  for  laser  illumination  and  atom 
introduction  into  the  trap. 

Overcoming,  or  at  least  reducing,  the  induced  kick 
is  the  easiest  part  of  this  problem.  By  using  a  detuning 
for  maximum  reflection,  the  kick  is  reduced  by  an  order  of 
magnitude  from  its  value  when  detuning  for  a  moderate  kick 
is  used.  Three  other  factors  are  introduced  by  using  the 
higher  detuning.  First,  the  maximum  normal  velocity  is 
raised,  deepening  the  trap.  Second,  the  length  of  time  of 
the  interaction  is  reduced  by  an  order  of  magnitude  due  to 
the  increased  potential.  This  reduces  the  chances  of  an 
unwanted  atomic  transition  out  of  resonance,  and  thereby 
improves  the  overall  efficiency  of  the  reflection.  Finally, 
the  gaussian  spread  of  velocities  is  reduced  by  the  use  of 
a  higher  detuning.  This  follows  from  the  Quantum  Considera¬ 
tions  section  of  the  Basic  Mechanism  discussion. 

Design  of  the  x  variation  of  the  trap  is  more  diffi¬ 
cult.  The  obvious  solution  is  to  extend  the  lateral  shape 
to  the  longitudinal  direction.  This  will  result  in  an 
elliptic  cup  which  would  return  all  atoms  to  their  original 
apex  position.  The  problem  with  this  approach  originates 
in  the  inherent  kick.  To  preserve  the  lateral  confinement 
the  laser  radiation  must  be  directed  along  the  x  axis.  This 
produces  a  kick  in  the  x  direction.  For  an  anticipated  size 
of  the  trap  on  the  order  of  h  =  10cm  ,  the  normal  velocity 
due  to  gravity  is  140cm/s  on  axis.  This  produces  a  kick  of 


5.2cm/s  on  each  bounce.  This  is  larger  than  the  sigma  for 
the  velocity  spread  (1.8cm/s),  but  the  problem  is  really  in 
the  consistency  of  direction  for  multiple  bounces.  The 
velocity  increments  due  to  the  kick  are  additive,  leading 
to  an  acceleration  out  of  the  trap,  whereas  the  gaussian 
spreading  results  in  no  net  acceleration. 

Two  alternative  approaches  present  themselves  for 
solution  of  this  net  acceleration.  The  first  is  to  switch 
the  direction  of  the  laser  on  a  time  interval  which  will 
catch  the  atoms  on  alternate  bounces.  This  would  be  tricky 
and  wasteful  of  laser  resources.  In  addition,  the  atoms 
have  a  different  time  between  bounces  dependent  on  their 
velocity  perpendicular  to  the  laser  axis.  Finally,  this 
approach  has  the  drawback  of  losing  atoms  which  strike  the 
surface  during  the  laser  direction  switch.  This  will 
quickly  lead  to  depletion  of  the  trapped  population  and, 
more  importantly,  damage  to  the  surface  and  consequent 
degradation  of  the  reflection  and  trapping.  The  second 
approach  is  to  tilt  the  x  variation  of  the  elliptical  cup  to 
allow  gravity  to  counteract  the  inherent  acceleration.  This 
approach  will  be  examined  in  detail. 

The  addition  of  a  tilt  to  the  ellipse  for  longitudi¬ 
nal  trapping  involves  the  addition  of  an  extra  term  in  the 
differential  equation  defining  the  surface.  The  equation 


becomes 


Unfortunately,  this  results  in  a  nonlinear  equation  which 
defies  analytic  solution.  Because  R  is  a  small  parameter  in 
the  proposed  application,  perturbation  theory  was  used  to 
obtain  an  approximate  form  for  the  longitudinal  shape.  The 
equation  can  be  written  (to  an  error  on  the  order  of 
R  x  2/4h2)  as 


z  =  Rx  +  £■  ♦  24  ♦  <8r2+P*4 
4h  4h2  32h3 


(34) 


The  purpose  of  the  shape  change  is  to  absorb  the 
additional  momentum  of  the  kick  and  to  redirect  it  in  the 
z  direction.  This  results  in  an  incremental  raising  of  the 
apex  position  with  no  change  in  the  longitudinal  position. 

The  effect  of  the  gaussian  spread  is  to  shift  the 
position  of  the  trajectory  apex.  This  causes  a  misalignment 
of  the  trajectory  with  the  design  parameters  of  the  trap. 

This  will  not  normally  lead  to  escape  from  the  trap,  but 
results  in  a  spread  in  the  population  available  for  observa¬ 
tion  at  the  apex  position.  The  position  spread  depends  in 
a  complicated  manner  on  the  initial  velocity  components  of 
the  atom,  and  how  long  it  stays  in  the  trap.  The  effect  of 
a  misalignment  tends  to  be  carried  through  subsequent  bounces 
and  is  slightly  enhanced  rather  than  mitigated.  Repeated 


simulations  have  shown,  however,  that  the  variation  is  on 
the  order  of  1cm  for  the  trap  parameters  given  after  ten 
bounces,  which  takes  2.5sec. 

The  reliable  position  of  the  atomic  apex  on  time 
scales  measured  in  seconds  is  at  the  crux  of  this  trap 
design.  For  a  trap  scale  of  h  =  10cm  ,  the  atoms  will  be 
at  the  apex  at  -0.27  second  intervals.  They  will  fall  1cm 
in  0.045  seconds.  Therefore,  fully  one  third  of  the  time 
an  atom  is  trapped,  it  is  in  the  top  10%  of  the  trap  with 
vertical  velocity  below  45cm/s.  Equipment  can  then  be 
trained  on  the  apex  position  to  use  slow  atoms  in  free 
flight,  without  laser  interactions  and  with  a  very  small 
spread  of  velocities. 

Laser  illumination  of  the  trap  is  somewhat  tricky. 

At  the  interior  interface,  the  radiation  must  be  aligned  in 
the  x  direction.  This  is  topologically  impossible.  The 
nearest  realizable  alternative  is  to  reflect  the  radiation 
along  the  interior  of  the  trap  in  a  meridian  flow.  That  is 
to  illuminate  the  trap  at  one  end  (y  =  0)  with  an  expanding 
wavefront.  Ignoring  the  tilt  correction,  if  the  focus  of 
the  laser  is  at  (/2h,0,h)  the  radiation  will  follow  the 
curve  of  the  bowl  around  to  (-/2h,0,h).  These  paths  result 
from  using  a  constant  thickness  of  dielectric  for  the  trap 
and  ignore  the  gaussian  propagation  anomalies  of  the  laser 
beam.  This  oversight  can  be  mitigated  by  cutting  a  small 
semi-circle  out  of  the  dielectric  around  the  laser  focus  and 


beveling  this  edge  for  normal  incidence.  The  correction  to 
this  procedure  necessitated  by  the  longitudinal  tilt  is 
straightforward.  The  radiation  will  still  follow  the  merid¬ 
ians  of  the  tilted  bowl. 

There  are  two  further  problems  with  the  above 

arrangement.  The  first  is  the  requirement  for  large  laser 

intensities.  At  the  center  ellipse  of  the  trap  (x  =  0)  the 

length  of  the  ellipse  is  approximately  hTr/3/2  .  For 

h  =  10cm  ,  this  gives  ~38.5cm.  For  dielectric  thickness  d, 

the  laser  covers  a  longitudinal  length  2d  on  each  reflection. 

Therefore,  the  laser  illumination  is  spread  over  an  area 
o 

~77d  cm  .  With  a  thickness  of  2mm  and  laser  intensity  of 
lwatt ,  a  sufficient  electric  field  can  be  generated.  Smaller 
values  of  d  will  result  in  unwanted  waveguide  effects.  The 
second  problem  involves  the  placement  of  observation  equip¬ 
ment.  This  can  be  accommodated  by  restricting  the  size  of 
the  trap  lateral  to  the  laser  propagation.  By  removing  two 
sections  from  the  top  of  the  trap  (10  degrees  meridionally), 
sufficient  room  is  provided  for  horizontal  observation  of 
the  apex  region  without  significant  loss  to  the  trapping 
potential  of  the  device.  This  has  the  advantage  of  reducing 
the  laser  power  required,  though  only  slightly. 

The  introduction  of  the  atoms  into  the  trap  is 
another  problem  which  should  be  addressed.  The  atoms  can 
enter  the  trap  through  a  small  hole  (~lmm),  low  on  the  side 
of  the  trap.  If  they  have  approximately  the  right  velocity 
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and  are  sufficiently  collimated  they  will  then  bounce 
freely  in  the  trap  until  they  spontaneously  shift  out  of 
resonance  or  pick  up  sufficient  horizontal  velocity  to 
escape.  Simulation  has  shown  that  the  spread  in  apex  posi¬ 
tion  makes  the  atoms  unsuitable  for  observation  after  10  to 
20  bounces.  This  should  be  before  either  of  the  other 
effects  removes  them  from  the  trap.  Each  atom  has  an  effec¬ 
tive  observation  time  on  the  order  of  one  second  spread  over 
a  lifetime  of  three  seconds. 

The  requirements  on  the  collimation  and  velocity  of 
the  injected  atoms  are  stringent.  However,  the  beam  nozzle 
of  the  preceding  section  and  deceleration  by  a  suitable  beam 
pipe  are  ideally  suited  to  provide  a  large  number  of  suitably 
prepared  atoms. 


III.  Analysis  and  Conclusions 

General  Considerations 

An  analysis  of  evanescent  resonant  reflection  and 
three  possible  applications  has  been  presented  in  the  pre¬ 
ceding  sections.  The  analysis  has  focused  on  the  production 
and  manipulation  of  very  well  collimated  and  very  slow 
(~100cm/s)  beams  of  neutral  particles.  The  questions  which 
should  be  asked  of  the  applications  presented  are  as  follows 
First,  are  the  proposed  applications  feasible  using  existing 
technology?  Secondly,  are  the  applications  useful  in  the 
study  of  neutral  particle  phenomena?  And  finally,  are  these 
applications  more  effective  than  existing  and  proposed 
methods  for  the  manipulation  of  atomic  beams?  These  are  the 
considerations  which  should  be  addressed  in  an  evaluation 
of  these  results. 

Criteria  for  Evaluation 

In  assessing  the  answers  to  the  three  considerations 
of  the  previous  section,  quantifiable  criteria  should  be 
adopted.  This  section  attempts  to  define  such  criteria. 

The  question  of  feasibility  involves  a  knowledge  of 
current  technology  in  tunable  laser  power.  In  addition,  the 
manufacture  of  precise  dielectric  surfaces  and  the  focusing 
of  the  available  laser  power  must  be  considered.  Finally, 
the  current  state  of  neutral  beam  production  should  be 


addressed  as  it  relates  to  the  input  parameters  for  any  of 
the  proposed  application.  Available  technology  is  the  best 
guarantor  of  feasibility. 

Current  tunable  laser  power  technology  is  available 
from  a  literature  search.  Powers  of  one  watt  are  common  in 
the  literature  of  sodium  experiments  (Refs  4-7).  One  of  the 
largest  power  levels  found  was  40kW  (pulsed)  for  a  tunable 
dye  laser  (Ref  32).  The  nominal  one  watt  power  level  used 
in  assessing  the  applications  in  this  thesis  is  easily 
attainable. 

The  ability  to  manufacture  precise  dielectric  sur¬ 
faces  does  not  appear  to  be  a  problem.  Current  lens  manu¬ 
facturing  techniques  can  be  used  to  form  the  necessary 
shapes.  Focusing  the  laser  power  onto  the  interface  is 
likewise  well  understood.  The  major  problems  are  waveguide 
and  diffraction  effects  which  are  combatted  by  using  thicker 
dielectric  slabs  or  more  laser  power. 

Neutral  beam  production  is  the  current  weak  link  in 

the  study  of  neutral  atom  characteristics.  Atomic  ovens 

produce  a  large  flux  of  atoms  into  a  2tt  solid  angle.  The 

problem  is  that  baffles  must  currently  be  used  to  effect 

collimation.  Beam  collimations  from  .01  to  .0005  radian  are 
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used  in  beam  experiments  (Ref  18).  These  equate  to  5x10 

_7 

to  1.25x10  of  the  atoms  issuing  from  the  oven.  This 
means  that  months  are  required  to  observe  the  equivalent  of 
minutes  of  oven  output . 
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The  usefulness  of  applications  must  be  based  on  the 
end  uses  for  the  atomic  beams.  The  major  uses  are  collision 
and  deflection  experiments  with  very  well  collimated  beams, 
and  atomic  spectroscopy  with  slow  beams  or  trapped  atoms. 

In  either  case,  the  criteria  reduce  to  the  spread  of  veloci¬ 
ties.  For  beam  collimation,  lateral  velocity  is  most 
important.  For  trapping,  all  components  of  velocity  are 
considered. 

The  effectiveness  of  beam  control  techniques  must 
be  judged  on  the  final  velocity  spread  and  the  required 
laser  power  and  oven  production.  This  is  the  final  crite¬ 
rion  for  beam  control  applications. 

Comparison  With  Existing  Results 

Deflection .  The  deflection  of  neutral  beams  was 
made  possible  only  recently  by  the  use  of  resonance-radiation 
pressure  (Refs  1-3;  18;  22).  The  spontaneous  force  was  used 
by  Scheider  et  al.  (Ref  18)  to  deflect  a  beam,  but  only 
marginal  success  was  achieved  (~28cm/s).  Cook  and  Bernhardt 
proposed  (Ref  22)  and  Arimondo  et  al .  confirmed  (Ref  3)  the 
possibility  of  using  the  induced  transition  rate  for  deflec¬ 
tion.  This  results  in  a  statistical  variation  with  ftt  as 
the  driving  variable.  Deflections  from  a  standing  wave 

pattern  using  this  scenario  can  be  very  high.  The  tangent 

2 

of  the  def1  ction  angle  varies  as  1/v  .  For  initial  veloci- 


Using  the  electromagnetic  mirror  of  this  thesis,  a 
1°  deflection  can  be  achieved  for  beam  velocities  below 
52000cm/s  (using  detuning  for  maximum  reflection  and  one 
watt  power).  Therefore,  the  evanescent  effect  is  20  times 
weaker  than  a  standing  wave  deflection.  The  advantage  of 
the  evanescent  reflection  is  its  sharply  controllable 
position  and  orientation. 

Focusing.  The  controllability  of  evanescent  reflec¬ 
tion  is  most  evident  in  the  focusing  of  a  diverging  atomic 
stream. 

The  current  method  of  focusing  a  neutral  atomic  beam 
is  to  co-propagate  a  gaussian  laser  beam  tuned  below  the 
resonant  transition.  This  causes  a  dipole  force  on  the 
atoms  which  deflects  them  toward  the  center  of  the  laser 
beam  (the  most  intense  position).  This  focusing  effect  has 
been  observed  by  Bjorkholm  et  al.  (Ref  5).  They  report  that 
the  atomic  trajectories  are  sinusoidal  at  a  frequency  of 
~4kHz  with  transverse  confinement  for  radial  velocities  below 
190cm/s. 

There  are  two  objections  to  this  method  of  focusing. 
The  first  is  that  the  beam  can  only  be  confined  along  the 
straight  path  of  the  laser.  For  most  applications  this  is 
only  a  minor  hindrance.  The  second  objection  is  that  the 
sinusoidal  trajectories  provide  good  containment  but  limited 
focusing  or  collimation.  Ignoring  the  gaussian  velocity 
spread,  the  radial  velocity  of  a  contained  atom  is  less  than 
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10%  of  its  maximum  value  for  6%  of  its  trajectory.  The 
radial  velocity  is  within  10%  of  its  maximum  value  for  29% 
of  the  time. 

The  enhanced  nozzle  of  this  thesis,  in  contrast, 
redirects  the  radial  components  of  a  divergent  atomic  beam. 
The  effect,  again  ignoring  the  gaussian  velocity  spread,  is 
to  produce  a  perfectly  collimated  beam.  For  other  applica¬ 
tions,  the  paraboloid  can  be  reshaped  into  an  ellipsoid. 
This  will  focus  a  large  portion  of  the  beam  onto  a  small 
spot.  No  comparable  effect  has  been  reported  in  the  study 
of  vacuum  resonant  radiation  pressure. 

Trapping.  The  current  literature  contains  many 
proposals  for  atomic  traps.  The  most  promising  design  is 
that  of  A.  Ashkin,  first  outlined  in  1978  (Ref  10).  This 
design  appears  capable  of  trapping  and  cooling  atoms  whose 
initial  velocity  components  are  lesfe  than  500cm/s  radially 

4 

and  2x10  cm/s  axially.  An  important  measure  of  trap  per¬ 
formance  is  the  Boltzmann  factor:  the  ratio  of  confining 
potential  well  depth  to  the  minimum  energy  attainable. 
Ashkin  calculates  a  Boltzmann  factor  of  ~100  for  his  trap. 
This  compares  favorably  to  the  Boltzmann  factor  of  ~1  for 
other  trapping  proposals. 

The  trap  described  in  this  thesis  has  a  Boltzmann 


factor  of  -2-10  depending  on  the  design  size  and  the  intro¬ 
duction  conditions  placed  on  the  atoms.  The  major  advantage 
of  this  trap  over  all  others  in  the  literature  is  the 


illumination-free  observation  of  the  atoms.  The  trap  has  a 
slightly  longer  containment  time,  a  ten  thousand  fold 
increase  in  the  laser-free  observation  time,  and  more  stable 
atomic  velocities  than  the  Ashkin  trap.  When  used  in  con¬ 
junction  with  an  enhanced  nozzle  and  beam  pipe  deceleration, 
either  trap  can  enjoy  a  greatly  enhanced  trapping  rate. 

Conclusions 

The  phenomena  of  evanescent  resonant  reflection  is 
worthy  of  further  theoretical  study  and  experimental  verifi¬ 
cation.  The  unique  properties  of  this  phenomena  show  great 
initial  promise  in  the  field  of  neutral  beam  control.  There 
are  distinct  possibilities  in  the  areas  of  controlled  deflec 
tion,  atomic  trapping,  and  particularly  in  the  collimation 
and  focusing  of  divergent  beams. 

Further  theoretical  study  should  concentrate  on  the 
strength  of  induced  diffusion  normal  to  the  reflecting 
surface  and  the  identification  of  additional  promising 
geometries  and  applications.  The  induced  diffusion  is  still 
the  major  limiting  factor  on  the  applicability  of  this 
technique.  For  the  parameters  of  this  thesis,  the  induced 
diffusion  precludes  specular  atomic  reflection  for  normal 
velocities  above  ~400cm/s. 

Experimental  verification  of  the  essential  features 
of  the  effect  is  also  needed.  The  magnitude  of  the  forces 


and  diffusion  coefficients  involved  can  be  verified  by 
conducting  fairly  simple  deflection  experiments.  The  con 
struct ion  and  testing  of  a  beam  pipe  should  also  yield 
results  which  would  confirm  (or  deny)  the  analysis  pre¬ 
sented. 
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Appendix  A 

Optics  Communication  Article 


The  genesis  of  this  thesis  was  an  idea  by  Dr.  Richard 
Cook  of  the  AFIT  staff.  This  idea  is  expressed  and  expounded 
in  the  following  article  accepted  for  publication  in  Optics 
Communications .  The  article  outlines  the  use  of  evanescent 
waves  for  the  manipulation  of  beams  of  neutral  atoms.  The 
quantum  mechanical  basis  for  the  force  laws  is  reviewed,  and 
conjectures  are  made  as  to  possible  applications.  The  author 
of  this  thesis  is  included  as  a  co-author  of  the  article 
because  of  efforts  in  verifying  the  derivation  and  assessing 
the  implications  for  neutral  beam  control. 

The  article  is  presented  as  written,  with  only 
necessary  modifications  to  the  form,  bibliography,  and 
citation.  The  salient  features  of  the  article  are  also 
presented  in  the  Basic  Mechanism  portion  of  section  II  of 
this  thesis.  Note  the  citations  are  to  the  bibliography 
on  page  52  while  the  equations  are  numbered  as  they  appear 
in  the  journal  article. 


AN  ELECTROMAGNETIC  MIRROR 
FOR  NEUTRAL  ATOMS 

Richard  J.  Cook 
Richard  K.  Hill 


Abstract 


When  light  is  totally  reflected  internally  at  a 
vacuum-dielectric  interface,  an  atom  in  the  thin  transmitted 
evanescent  wave  experiences  a  radiation  force.  For  light 
tuned  above  the  transition  frequency  of  a  two-level  atom, 
this  force  tends  to  repel  the  atom  from  the  dielectric 
surface;  and  hence  the  internally  illuminated  surface  acts 
as  a  mirror  for  slow  neutral  atoms.  This  paper  presents 
the  first  analysis  of  this  atomic  reflection  process  and 
suggests  that  the  effect  can  be  used  to  trap  slow  atoms  or 


to  focus  a  slow  atomic  beam. 
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In  recent  years,  a  number  of  techniques  have  been 
proposed  for  controlling  the  motion  of  atoms  or  molecules  by 
means  of  the  radiation  force  of  laser  light  propagating 
freely  in  vacuum  [5;  9-13;  17;  20-23;  27;  28].  This 
paper  describes  a  new  scheme  for  controlling  atomic  motion 
which  makes  use  of  the  radiation  force  exerted  by  the  thin 
evanescent  wave  that  is  generated  on  the  surface  of  a 
dielectric  medium  when  laser  light  is  totally  reflected 
internally  at  that  surface.  Specifically,  we  show  that, 
when  the  light  is  tuned  above  an  atomic  resonant  frequency, 
a  slow  atom  is  reflected  by  the  evanescent  wave  without 
making  contact  with  the  material  surface. 

If  a  plane  electromagnetic  wave  of  frequency  oj 
propagates  in  a  medium  of  refractive  index  n  and  is  totally 
reflected  internally  as  depicted  in  Figure  8,  then  the 
electric  field  of  the  transmitted  evanescent  wave  takes  the 
form  [29] 

2(x,t)  =  ££e-aycos(wt-kx)  Cl) 

where  §  is  a  unit  polarization  vector, 

a  =  w(n2sin29  -  l)*/c,  (2) 

k  =  cunsin8/c  (3) 

0  is  the  angle  of  incidence  of  the  wave,  and  £  is  the  wave 
amplitude  at  y  =  0  .  The  condition  for  total  internal 
reflection  is  0  <  0  =  sin“^(l/n)  .  The  evanescent  wave 

C 
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propagates  parallel  to  the  dielectric  surface  and  is  sub¬ 
stantial  only  within  a  few  wavelengths  of  the  surface, 


except  for  8  *  8C  . 

Now  an  atom  in  the  evanescent  wave  experiences  a 
radiation  force  due  to  momentum  transfer  from  the  wave.  For 
a  two-level  atom  with  transition  frequency  cuo  and  a  near¬ 
resonant  field  (w  »  0)o),  the  radiation  force  has  a  component 
[30] 


2afi(A-kvx)02(y) 

4(A-kvx)2+A2+2fi2(y) 


(4) 


normal  to  the  dielectric  surface  associated  with  the  gradient 
of  the  field  amplitude  (the  dipole  force),  and  a  component 


F 


x 


hk AQ2 ( y ) 

4 ( A-k v  ) 2+ A2+202 ( y ) 


(5) 


parallel  to  the  surface  (the  spontaneous  force),  where 
A  =  to— to  is  the  detuning  frequency,  v  is  the  x-component 

O  X 

2  3 

of  the  atomic  velocity,  A  =  4y  u>0/3hc  is  the  Einstein 
spontaneous  emission  coefficient,  and  ft(y)  =  y£exp(-ay)/li 
is  the  on-resonance  Rabi  frequency  of  the  two-level  atom 
(u  is  the  dipole  transition  moment).  Equations  (4)  and  (5) 
are  valid  when  the  Rabi  frequency  0(y)  and  the  Doppler  shift 
kv  are  slowly  varying  functions  of  time  on  the  time  scale 
of  the  atomic  relaxation  time  x  =  1/A  [30].  For  simplicity, 
the  following  discussion  is  limited  to  the  case  in  which  the 


detuning  A  is  positive  (w  >  u>Q )  and  is  large  compared  to  the 
Doppler  shift  kvx.  Then,  to  a  good  approximation,  Eqs  (4) 
and  (5)  reduce  to 


p  =  2afiAfl2(y ) 

y  4A2+A2+202( y )  (6) 

F  =  RkAfl2(y ) 

X  4A2+A2+2fl2(y)  (7) 


Because  A>0  ,  F  tends  to  repel  the  atom  from  the  dielectric 

y 

surface.  The  normal  force  F  ,  taken  alone,  is  derivable 
from  the  potential  energy 


V(y)  =  (l/2)RAln[l+202(y)/4A2+A2)] 


(8) 


Therefore,  an  atom  of  mass  II  initially  moving  toward  the 

surface  with  velocity  vy  is  reflected  by  the  thin  evanescent 
o 

wave  if  JMv  <  V(0)  ,  i.e.,  the  maximum  reflected  velocity 

y 

is 

max  .  lt  2vV  )>*  (9) 

y  M  h2( 4A2+A2) 


Next  we  derive  the  law  of  reflection  for  the  atomic 
trajectory.  From  Eqs  (6)  and  (7)  the  ratio  of  the  tangential 
to  the  normal  force  is  seen  to  be 


R  » 


kA 
2ct  A 


nAsinQ 

2A(n2sin20-l)* 


(10) 
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This  ratio  is  independent  of  position.  Therefore,  the  in¬ 
crements  of  velocity,  Av  and  Av  ,  experienced  by  the  atom 

x  y 

during  the  reflection  process  must  have  the  same  ratio 

(Av  =  RAv  ) ,  and  since  in  reflection  Av  =  -2v  ,  where 

x  y  y  y 

vy  is  the  initial  (negative)  normal  velocity,  we  have 

Av  =  -2Rv  .  Suppose  now  that  the  atom  moves  in  the  xy 
x  y 

plane  as  illustrated  in  the  Figure.  Then  from  the  geometry 
of  the  initial  and  final  velocity  vectors,  it  is  easy  to 
show  that  the  atomic  trajectory  obeys  the  law  of  reflection 

tan<J>r  =  tan<j>i  -  2R  (11) 

where  the  angle  of  incidence  4k  and  the  angle  of  reflection 
4>r  are  algebraic  quantities  measured  in  the  sense  indicated 
in  Figure  1.  Note  that  a  positive  4k  indicates  that  vx  is 
initially  opposite  to  the  propagation  direction  of  the 
evanescent  wave.  Equation  (11)  has  the  solution  4>r  =  -4k 
for  4^  ■  tan~*R  .  For  this  angle  of  incidence,  the  atom  is 
retroreflected.  When  R  is  negligibly  small  (R<<1),  Eq  (11) 
reduces  to  the  simple  reflection  la*"  4>r  =  4^  •  According 
to  Eq  (10),  this  occurs  when  the  detuning  A  is  much  larger 
than  the  Einstein  A-coef  f  icient ,  except  for  8  **  0C  = 
sin-1(l/n)  .  It  should  be  noted  that  if  the  atom  has  a 
component  of  velocity  in  the  z-direction,  the  planes  of  in¬ 
cidence  and  reflection  of  the  atomic  trajectory  are  not, 
in  general,  the  same  because  of  the  impulse  delivered  to 
the  atom  in  the  x-direction.  However,  as  R+0  these  planes 


become  coincident  and  the  simple  reflection  law  4>r  =  <J>^ 


is  again  obtained. 


To  estimate  the  atomic  velocity  v^  that  can  be 
reflected  in  practice,  we  suppose  that  a  z-polarized  elec¬ 
tromagnetic  wave  of  intensity  I  enters  the  dielectric  norm¬ 
ally  through  a  plane  surface  not  shown  in  Figure  8.  By 
a  straightforward  application  of  Fresnel's  reflection 
formulas  [29]  ,  it  is  found  that  the  square  of  the  evanescent 
wave  amplitude  is  £2  =  1287rn2cos20I/(n2-l)(n+l)2/c  .  Let 
the  intensity  and  angle  of  incidence  of  the  wave  by  1=1 
watt/cm2  and  9  =  60°  ,  respectively,  and  let  the  refractive 
index  of  the  dielectric  be  n  =  1.5  .  Then  for  reasonable 

atomic  parameters  (M  =  4x10  g  ,  wQ  =  3x10  s'  ,  A  =  10°s 

8  —1 

and  detuning  A  =  5x10  s"  ),  one  obtains  from  Eq  (9)  the 
max 


In  a  strong  field  (0Q  =  y£/h>>A) t 


value  =  llOcm/s  . 

vmax  assumes  maximum  value  v™3*  *  0.75(hOQ/M)^  for 

detuning  A  s  0.56flQ  ,  and  since  a  1^  ,  this  implies 

nisix  ^ 

that  the  upper  limit  on  v  is  proportional  to  I4,  which 
is  a  rather  slowly  increasing  function  of  intensity. 
Therefore,  the  available  tunable  CW  laser  power  (~lwatt) 

limits  v^3-*  to  values  of  order  100-400cm/sec  for  mirror 

2  2  max 

areas  ranging  from  1cm  to  1mm  .  However,  values  of  vy 

4 

approaching  typical  thermal  velocities  (-10  cm/s)  can  be 
achieved  for  very  short  time  intervals  by  use  of  pulsed 
laser  radiation. 
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The  above  considerations  suggest  a  new  type  of  trap 

for  slow  neutral  atoms.  Suppose  that  laser  radiation  is 

totally  reflected  by  each  of  the  six  faces  of  a  cubical 

cavity  in  a  dielectric  medium.  Then  a  slow  atom  inside  the 

cavity  will  be  reflected  by  the  evanescent  waves  on  the 

cavity  walls  without  making  thermal  contact  with  the 

dielectric  itself.  In  this  way  a  cold  atomic  vapor 
-4 

(T  <»  5x10  K)  might  be  held  in  the  cavity  for  a  time 
sufficiently  long  to  perform  essentially  Doppler-free 
spectroscopy.  An  advantage  of  this  type  of  trap  is  that  the 
atoms  interact  with  the  radiation  only  during  the  brief 
reflection  processes,  and  consequently  the  quantum  fluctua¬ 
tions  of  the  radiation  force  [14;  31],  which  tend  to  heat 
the  atomic  vapor,  are  much  less  of  a  problem  than  in  those 
schemes  in  which  the  atoms  are  continuously  illuminated. 

Other  possible  applications  of  the  reflection  pro¬ 
cess  immediately  come  to  mind.  For  example,  it  is  evident 
that  the  dielectric  surface  in  Figure  8  can  be  given  a  small 
curvature  to  produce  a  concave  mirror  for  neutral  atoms. 

Such  a  mirror  could  be  used  to  focus  an  atomic  beam  of  slow 
atoms  or  to  image  atoms  that  are  scattered  by  a  small  object. 
Moreover,  a  complete  "optics"  of  atomic  beams  could  be  con¬ 
structed  on  the  basis  of  such  reflecting  elements.  But  in 
this  connection,  it  is  very  important  to  note  that  gravity 
will  play  an  important  role,  since,  for  the  velocities  of 


interest  ( -100cm/ sec ) ,  the  atomic  trajectories  are  bent 
significantly  by  gravity  over  a  distance  of  a  few  centi¬ 
meters. 
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Appendix  B 
BASIC  Program 


General 

In  an  effort  to  understand  and  examine  the  basic 
mechanism  of  resonant  reflection,  and  to  test  possible 
applications,  the  following  computer  program  was  written. 

The  language  used  is  BASIC  as  implemented  on  an  H89  using 
CP/M.  The  program  was  made  as  flexible  as  possible  and 
incorporates  several  user-defined  routines. 

The  user  can  define  any  mirror  geometry  within 
reason.  As  listed,  the  program  defines  the  mirror  surface 
ZSURF  as  a  function  of  X  and  Y.  Complete  three  dimensional 
surfaces  are  possible  with  slight  modifications,  as  was 
used  to  test  the  Enhanced  Nozzle  application.  This  simply 
involves  the  use  of  in  and  out,  instead  of  above  and  below, 
in  the  definition  of  the  surface  and  the  intercept  location 
routine.  For  the  sake  of  execution  speed,  a  separate  routine 
calculates  the  mirror  normal  for  orientation  of  the  reflec¬ 
tion. 

Three  dimensional  laser  variations  are  also  definable 
in  the  WAVEFRONT  routine.  The  subroutine  calls  are  arranged 
so  that  the  orientation  of  the  surface  is  available  when  the 
wavefront  routine  is  called.  This  is  extremely  useful  in 
obtaining  standard  reflection  parameters  to  investigate  a 


surface  geometry  without  concurrent  solution  of  the  full 
laser  variation  required.  This  greatly  speeds  the  task  of 
investigating  myriad  and  unusual  mirror  shapes. 

The  BOUNCE  routine  consists  of  three  sections  which 
may  be  mixed  as  desired  by  inserting  jumps  over  unwanted 
sections.  The  first  section  calculates  a  normal  elastic 
bounce.  This  was  useful  in  testing  certain  notions  of 
mirror  properties  without  the  added  complication  of  the 
induced  kick.  The  second  section  calculates  a  bounce 
including  the  kick.  On  each  bounce  it  uses  the  current 
evanescent  wave  strength  and  direction  to  compute  the  bounce 
and  kick.  It  also  checks  that  the  normal  velocity  is 
within  allowable  parameters.  The  third  section  calculates 
a  random  variation  to  simulate  the  gaussian  spread  in 
velocities  due  to  the  spontaneous  nature  of  the  force. 

The  gaussian  spread  simulation  is  worthy  of  further 
note.  For  each  direction  in  the  interaction  reference 
frame,  a  random  number  is  chosen  between  zero  and  one.  This 
number  is  then  used  as  the  target  in  an  iterative  search 
which  uses  a  ninth  order  polynomial  expansion  of  the  error 
function.  The  result  is  a  random  variable  with  a  gaussian 
probability  distribution.  The  root-mean-square  deviation 
of  the  spread  is  then  calculated  based  on  the  normal  velocity 
and  evanescent  wave  parameters.  Multiplication  of  this  sigma 
by  the  calculated  random  numbers  gives  a  representative 
velocity  shift  for  each  direction.  The  velocity  increments 


are  added  in  the  interaction  frame,  but  the  effect  in  any 
frame  is  the  same  gaussian  spread  in  velocity. 

Program  Use 

When  the  program  is  run,  it  first  checks  the  status 
of  the  line  printer  and  asks  for  it  to  be  turned  on  if 
desired.  The  program  then  prints  out  messages  to  identify 
the  mirror  shape,  laser  incidence  parameters,  and  type  of 
bounce  the  program  is  currently  using.  Each  of  the  messages 
can  be  reprogrammed  to  aid  in  the  identification  of  results. 

The  user  is  then  asked  for  the  starting  height  of  the 
atom  and  the  initial  x,  y,  and  z  components  of  its  velocity. 
The  initial  values  of  x  and  y  are  currently  set  a  0  (on  axis), 
to  preclude  entry  each  time  the  program  is  run.  This  is 
easily  changed  to  incorporate  more  flexibility,  and  all  six 
components  of  the  atom's  position  and  velocity  are  printed  in 
any  event. 

The  next  step  in  the  program  is  the  calculation  of 
the  apex  position  and  energy.  This  serves  as  a  reference 
and  is  printed  after  each  bounce.  The  program  then  proceeds 
to  find  the  next  bounce  position.  It  prints  the  elapsed 
time  since  the  initial  conditions,  the  bounce  position,  and 
the  number  of  the  bounce. 


The  program  then  prints  the  interaction  parameters 
at  the  position  of  the  bounce.  These  are  the  LASER  inten¬ 
sity  in  watts,  WX  WY  and  WZ:  the  normalized  components 


of  the  incident  wave  vector,  THETA:  the  laser's  angle  of 
incidence,  R:  the  kick  ratio,  and  IP  and  IC:  the  normalized 
components  of  the  evanescent  wave.  These  are  defined 
Parallel  and  Complimentary  to  the  projection  of  the  atomic 
trajectory  on  the  surface-.  Finally  the  DETUNING,  including 
the  Doppler  contribution,  is  printed  in  units  of  Mhz. 

The  next  output  consists  of  SIGMA:  the  gaussian 
spread  parameter,  DN  DP  and  DC:  the  random  velocity 
deviations  Normal  to  the  surface  and  Parallel  and  Compli¬ 
mentary  to  the  trajectory  projection.  The  normal  velocity 
of  the  atom  is  then  printed  for  comparison  with  the  maximum 
normal  velocity  permitted  by  the  interaction  parameters. 

The  SPEED,  and  the  X  Y  and  Z  components  of  velocity 
before  and  after  the  bounce  are  then  output.  This  allows  a 
direct  evaluation  of  the  effects  of  the  interaction. 

Finally,  the  apex  energy  and  position  are  printed  to  show 
any  shifts. 
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10  REM  TRAP 

20  PI-3.141592653000001# 

30  0=29979250000#  s  REM  C6S  UNITS 

40  GRAV— 980 . 665 
50  HBAR-1 .0545*10A-27 

60  OMEGAO— 5 . 08985E+ 1 4  s  REM  SODIUM  3S1/2  TO  3P3/2 
70  MASS— 3 . 81 729E-23  s  REM  SODIUM 

80  RANDOMI ZE  23382  :  REM  INITIALIZE  FOR  GAUSSIAN  SPREAD 

90  BOUNCE-O 

100  FORM 1 *="*##.###  &#####.##  &#####.## 

&#####.##  &##" 

110  F0RM2*-"  &#####.##  &#####.##  &#####.##  &#####.## 

&####.##" 

1 20  F0RM3*- *&##### . ##&##### . ##&##### . ##" 

130  F0RM4*- "&#####.##  &####.##" 

140  PR=INP<230)  s  REM  CHECK  PRINTER  STATUS 

150  PR-PR16 

160  PRINTERS PR/2-PR2) #2 

170  IF  PRINTERS  OR  FLAGPR=1  THEN  230 

180  PRINT  "Please  Turn  on  the  Printer  <if  desired) ,  and 
press  RETURN" 

190  A*= INKEY* 

200  IF  LEN< A*)=0  THEN  190 
210  FLAGPR-1 

220  GOTO  140  :REM  END  OF  PRINTER  CHECK 

230  TIME-0 

240  XO-O 

250  Y0—0 

260  PRINT 

270  IF  NOT  PRINTER-1  THEN  340  *  REM  NO  PRINTER 
280  GOSUB  1430  :  REM  SURFACE  ROUTINE 
290  LPRINT  ZSURF* 

300  GOSUB  1760  s  REM  WAVEFRONT  ROUTINE 
310  LPRINT  WAVE* 

320  GOSUB  2580  s  REM  BOUNCE  ROUTINE 
330  LPRINT  BOUNCE* 

340  INPUT | "  ZOs  ",Z0 
350  INPUT; "  VX:  " ,VX 
360  INPUT;"  VYs  " ,VY 
370  INPUT  "  VZs  ■ ,VZ 

380  IF  NOT  PRINTER-1  THEN  430  :  REM  NO  PRINTER 
390  LPRINT 

400  LPRINT  USING  F0RM3* ; "  XO *  ";X0;"  YO :  ";Y0;"  ZOs  ";Z0 
410  LPRINT  USING  F0RM3*;"  VXs  ";VX;"  VYs  ";VY;"  VZ s  ■ ;VZ 
420  LPRINT 

430  GOSUB  1340  s  REM  TRAJ  ROUTINE 
440  BSUBX-VX 
450  BSUBY-VY 
460  BSUBZ— VZ 

470  GOSUB  2120  s  REM  APEX  ROUTINE 
480  T— 0 


490  DT-.01 

500  GOSUB  1340  i  REM  TRAJ  ROUTINE 
510  GOSUB  1430  :  REM  MIRROR  SURFACE 
520  DZ-Z-ZSURF 
530  T*T+DT 
540  DZOLD-DZ 

550  GOSUB  1340  :  REM  TRAJ  ROUTINE 

560  GOSUB  1430  t  REM  MIRROR  SURFACE 

570  DZ=Z-ZSURF 

580  IF  DZ> .01  THEN  530 

590  IF  ABS<DZ) < .01  THEN  670 

600  DT=-DT*DZ/<DZ-DZOLD> 

610  T*T+DT 
620  DZOLD=DZ 

630  GOSUB  1340  :  REM  TRAJ  ROUTINE 
640  GOSUB  1430  :  REM  MIRROR  SURFACE 
650  DZ*Z-ZSURF 

660  GOTO  590  s  REM  NEXT  STEP  IN  INTERCEPT  SEARCH 
670  TIME=TIME+T 
680  BOUNCE=BOUNCE+l 
690  PRINT 

700  PRINT  USING  FORM1*  j  •TIME®*"  ;TIME;  *  X=';X;"  Y=*;Y;"  Z**  ; 

Z; ‘BOUNCE  *  ■ ;  BOUNCE 

710  IF  NOT  PRINTER=1  THEN  740  s  REM  NO  PRINTER 
720  LPRINT 

730  LPRINT  USING  FORM1*  ,■  "TIME*"  ;TIME  5  “  X=';X;-  Y=-;Y;-  Z*' 

Zj -BOUNCE  »  •; BOUNCE 
740  GOSUB  830  :  REM  BOUNCE  ROUTINE 

750  GOSUB  2120  :  REM  APEX  CALCULATION 
760  XO=X 
770  Y0=Y 
780  Z0=Z 
790  VX=BSUBX 
800  VY»BSUBY 
810  VZ=BSUBZ 

820  GOTO  480  s  REM  RETURN  TO  FIND  NEXT  INTERCEPT 
830  REM  BOUNCE 

840  GOSUB  1580  :  REM  MIRROR  NORMAL 
850  GOSUB  2040  s  REM  TRAJ  VECTOR 
860  GOSUB  1760  s  REM  WAVEFRONT  ROUTINE 

870  PRINT  USING  FORM 2* ; ‘ LASER®  • jINTEj'WX** ;WSUBXj"WY*- ; 
WSUBYj-WZ-' ;WSUBZ 

880  IF  NOT  PRINTERS  THEN  900  :  REM  NO  PRINTER 
890  LPRINT  USING  FORM 2* ; ‘LASER-  • 5 INTE ; "WX-* }WSUBX; "WY-- j 
WSUBYj'WZ-' jWSUBZ 

900  REM  TRANSFORMATION  TO  N,C,P  COORDINATES 
910  ADOTN=NORMX#WX+NORMY«NVY+NORMZ«WZ 
920  ASUBX-NVX 
930  ASUBY=N/Y 
940  ASUBZ"NVZ 

950  REM  DETERMINE  OF  P(ARALLEL)  fiHD  C<OMPLIMENTARY)  VECTORS 
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960  THETAI=1 

970  IF  NOT  (ADOTN  >-l)  THEN  THETAI-0 

980  IF  THETAI  =  0  THEN  1000  s  REM  AVOID  NEGATIVE  SQUARE  ROOT 
990  THETAI=ATN<-SQR<1-AD0TNA2>/AD0TN> 

1000  ASUBN— SPEED*COS<THETAI> 

1010  ASUBP=SPEED*SIN<THETAI > 

1020  ASUBC-0 

1030  IF  NOT  < THETAI  =  0)  THEN  1110  :  REM  NORMAL  INCIDENCE 
1040  CSUBX-0  s  CSUBY-1  s  CSUBZ-0 
1050  N0RMPERP=SQR<N0RMYA2+N0RMZA2> 

1060  IF  NORMPERP— 0  THEN  1140  *  REM  NORM  IS  THE  X  AXIS 
1070  CSUBX-0 

1080  CSUBY-NORMZ/NORMPERP 
1090  CSUBZ=-NORMY/NORMPERP 
1100  GOTO  1140 

1110  CSUBX»<ASUBY*NORMZ-ASUBZ*NORMY)/SIN<THETAI ) 

1120  CSUBY=<ASUBZ*NORMX-ASUBX*NORMZ> /SIN (THETAI) 

1 1 30  CSUBZ«<ASUBX*NORMY-ASUBY*NORMX>/SIN(THETAI  ) 

1140  PSUBX=NORMY*CSUBZ-NORMZ*CSUBY 

1150  PSUBY=NORMZ*CSUBX-NORMX*CSUBZ 

1160  PSUBZ=NORMX*CSUBY-NORMY*CSUBX 

1170  GOSUB  2260  :  REM  IMPULSE  DIRECTION  ROUTINE 

1180  GOSUB  2580  :  REM  BOUNCE  ROUTINE 

1190  REM  TRANSFORMATION  BACK  TO  X,Y,Z 
1 200  BSUBX=BSUBN*NORMX+BSUBC*CSUBX+BSUBP*PSUBX 
1210  BSUBY-BSUBN*NORMY+BSUBC*CSUBY+BSUBP*PSUBY 
1 220  BSUBZ=BSUBN*NORMZ+BSUBC*CSUBZ+BSUBP*PSUBZ 
1 230  SPEEDB-SQR  <  BSUBXA  2+  BSUBYA  2+ BSUBZ  A  2  > 

1240  A*='SPEEDA=* 

1250  BS=*  SPEEDB— * 

1260  C*=-VNORMMAX=  ■ 

1270  DS^'VNORM™ 

1280  PRINT  USING  F0RM2S {AS {SPEED {"VX—*  jVXj’VY*" ;VY;*VZ— * ; 
DZDT;CS;LNORMMAX 

1290  PRINT  USING  F0RM2S ; BS ; SPEEDB ; *  BX— ■ ; BSUBX ; * BY** ; BSUBY ; 
•BZ»' {BSUBZ 

1300  IF  NOT  PRINTER-!  THEN  1330  :  REM  NO  PRINTER 
1310  LPRINT  USING  F0RM2S ; AS { SPEED; BVX»* fVX; "VY-* ;VY{ "VZ-" ; 
DZDT ;CS;LNORMMAX 

1320  LPRINT  USING  F0RM2S ;BS{ SPEEDB {"BX—* {BSUBX {"BY—" {BSUBY; 

•BZ-- {BSUBZ 
1330  RETURN 

1340  REM  TRAJ  SUBROUTINE 
1350  REM  GIVEN  X0,Y0fZ0 
1360  REM  GIVEN  VX,VY,VZfT 
1370  REM  RETURNS  X,Y,Z,OZDT 
1380  X— X0+VX*T 
1390  Y-Y0*VY»T 
1400  Z— Z0+VZ#T-GRAV*T*T/2 
1410  DZDT— VZ-GRAV*T 
1420  RETURN 


1430  REM  MIRROR  SURFACE 
1440  REM  GIVEN  X,Y 
1450  REM  RETURNS  ZSURF 
1460  HE I 6HT= 1 0 
1470  PHI=ATN< .0187) 

1480  PHI D-PHI# 180/PI 
1490  Q-TAN(PHI) 

1500  ZSURF=X*Q 

1510  ZSURF=ZSURF+XA 2/4/HEIGHT 

1520  ZSURF=ZSURF+Q*XA3/4/HEIGHTA2 

1 530  ZSURF-ZSURF+  < 8*QA2+ 1 > *XA4/32/HEI GHTA3 

1560  ZSURF*— *Z  *  Four  Term  Expansion  of  Longitudinal  Trap 

1570  RETURN 

1580  REM  MIRROR  NORMAL 

1590  REM  GIVEN  X,Y 

1600  REM  RETURNS  NORMX ,N0RMY ,N0RMZ 

1610  G=TAN(PHI> 

1620  DZDX-Q 

1630  DZDX-D2DX+X/2/HE1GHT 
1640  DZDX=DZDX+3*Q*XA2/4/HEIGHTA2 
1650  DZ  DX=DZ  DX+  <  8*0 A  2+ 1 ) *XA  3/8/HE I GHTA 3 
1680  DZDY—0 

1690  GOTO  1710  t  REM  SKIP  LATERAL  FOCUS 

1700  DZDY=ALPHA*Y/2/ <  HEI GHT-Z  >  i  REM  LATERAL  FOCUS  2 
1710  L0RM=SQR(1+DZ0XA2+DZ0YA2> 

1720  NORMX— DZDX/LORM 
1730  NORMY=-DZDY/LORM 
1740  NORMZ=l/LORM 
1750  RETURN 

1760  REM  WAVEFRONT  ROUTINE 
1770  REM  GIVEN  X,Y,Z 

1780  REM  RETURNS  WSUBX,  WSUBY,  WSUBZ ,  WSTRENGTH ,  INDEX, 
DELTA,  EINA 
1790  THETAZ=ATN<DZDX> 

1800  THETAZ-THETAZ* 180/PI 
1810  THETAWAVE-45+THET AZ 
1820  THETAZ RAO— THETAWAVE*P 1/180 
1830  WSUBZ-S IN (THETAZ RAD) 

1840  WSUBX«C0S( THETAZ RAD) 

1850  WSUBY— 0 

1860  WAVE*-*45  Degree  Incidence  to  Surface  Upward* 

1870  GOTO  1940  s  REM  SKIP  PIPE  WAVEFRONT 

1880  THETAZ-45+PHIRAMPD 
1890  THETAZ RAD-THETAZ*PI/ 180 
1900  WSUBZ-C0S( THETAZ RAD) 

1910  WSUBX— SI N(THETAZRAD) 

1920  WSUBY-0 

1930  WAVE*- *WSUBZ-C0S( THETA)  WSUBX— SI N(THETA)  WSUBY-0 
THETA-' + ST R*< THETAZ > 

1940  INTENSITY-1 


l  REM  WATTS 


1950  INTENSITY*INTENSITY*10A7  :  REM  CGS  UNITS 
1940  INDEX=1 .5 

1970  DELTA=5*10A9  s  REM  DETUNING 

1980  OMEGA* 0MEGA0+ DELTA 

1990  EINA*4.22E+07  :  REM  SODIUM 

2000  THETAMIN=ATN< 1/INDEX/SQR< 1 -1/INDEXA2> > 

2010  THED=THETAM IN* 180/PI 
2020  INTE= INTENSITY/1 0A 7 
2030  RETURN 

2040  REM  TRAJ  VECTOR 

2050  REM  GIVEN  X,Y,Z,T 

2040  REM  RETURNS  SPEED,  NVX,NVY,WZ 

2070  SPEED*SQR<VXA2+VYA2+DZDTA2> 

2080  NVX=VX/SPEED 
2090  NVY=VY/SPEED 
2100  NVZ=DZDT/SPEED 
2110  RETURN 

2120  REM  APEX  CALCULATION 

2130  REM  GIVEN  B  VECTOR  COMPONENTS 

2140  REM  RETURNS  XAPEX , YAPEX , ZAPEX , ENERGYAP 

2150  TAPEX=BSUBZ/GRAV 

2140  T I MEAPEX=T I ME+TAPEX 

2170  XAPEX=X+BSUBX*TAPEX 

2180  YAPEX*Y  +  BSUBY *T APEX 

2190  ZAPEX=Z+GRAV*TAPEXA  2/2 

2200  SPEEDAPEX=SQR<  8SU8XA2+BSUBYA2) 

22 1 0  ENERGYAP*  <  SPEEDAPEXA 2/2+  GRAV*ZAPEX > / 1 0  0  0 
2220  PRINT  USING  F0RM2*; "ENERGY*- ; ENERGYAP ; ■  X*" ; XAPEX; 
•  Y*" j YAPEX j "  Z="  ; ZAPEX 

2230  IF  NOT  PRINTER*!  THEN  2250  :  REM  NO  PRINTER 
2240  LPRINT  USING  F0RM2*; "ENERGY* ENERGYAP ; "  X*" ; XAPEX 
"  Y*" j YAPEX; "  Z*" ; ZAPEX 
2250  RETURN 

2240  REM  IMPULSE  DIRECTION 

2270  REM  GIVEN  NORM,  CVECT,  PVECT,  WVECT 

2280  REM  RETURNS  IMPSUBP,  IMPSUBC,  RATIO 

2290  NDOTW*NORMX*WSUBX+NORMY*USUBY+NORMZ*WSUBZ 

2300  THETA-0 

2310  IF  NDOTW=0  THEN  2330  s  REM  NORMAL  LASER  INCIDENCE 
2320  THETA*ATN<SQR< 1-ND0TWA2)/ND0TW) 

2330  WPERX*WSU8X-N0RMX*C0S<THETA) 

2340  WPERY*«SUBY-NORMY*COS<THETA) 

2350  WPERZ*< JSUBZ -NORMZ  *COS <  THET A  > 

2340  IMPSUBP*<WPERX*PSUBX+WPERY*PSUBY+WPERZ*PSUBZ )/ 
SIN<THETA> 

2370  I MPSUBC* < WPERX*CSUBX*WPERY*CSUBY +WPERZ *CSUBZ > / 
SIN<THETA> 

2380  RATI 0-INDEX*EINA*S I N< THETA) /2/DELTA/ 

SQR( INDEXA2*SIN< THETA )A 2-1 > 

2390  THETD-THETA* 180/PI 


2400  REM  TEST  ON  EVANESCENT  WAVE 

2410  ALPHA2 OMEGA* SQR< INDEXA2*SIN<THETA> A2-l )/C 

2420  KUAVE=OMEGA*INDEX*SIN< THETA >/C 

2430  THETAMIN=ATN< 1/INDEX/SQR< 1 -1/INDEXA2>  > 

2440  ESQUARE=1 28*PI * INDEXA 2*C0S< THETA) A 2* INTENSITY/ 

< INDEXA2-1 >/< INDEX* 1 ) A2/C 
2450  MU=SQR< 3/ 4*EINA*HBAR> *SQR< < C/OMEGAO ) A3> 

2440  RABI=MU*SQR<ESQUARE>/HBAR 
2470  DDT2DELTA+KUAVE*ASUBP* IMPSUBP 
2480  TERM=2*MUA2/HBAR 
2490  TERM=TERM*ESQUARE/HBAR 
2500  TERM=TERM/<4*DDTA2+EINAA2> 

251 0  VNORf*tiAX=SQR<  HBAR*DDT/MASS*LOG<  1 +TERM)  > 

2520  PRINT  USING  F0RM2S; * THETA2  * ;THETD ; ■  R=* jRATIO; ■ IP2* ; 
IMPSUBP; 

2530  PRINT  USING  F0RM4S; " I C2‘ } I MPSUBC; 'DETUNING*  *;DDT/10A4 
2540  IF  NOT  PRINTER2!  THEN  2570  !  REM  NO  PRINTER 
2550  LPRINT  USING  F0RM2S ; "THETA2  ■ ;THETD; "  R2* ;RATIO; ■ IP2" ; 
IMPSUBP; 

2540  LPRINT  USING  F0RM4S; * IC2* ; I MPSUBC; “DETUNING2  "; DDT/1 0A4 
2570  RETURN 

2580  REM  BOUNCE  ROUTINE 

2590  BOUNCES* Bounce  without  Impulse* 

2400  BSUBN2 -ASUBN 
2410  BSUBP=ASUBP 
2420  BSUBC=ASUBC 

2430  BOUNCES2* Bounce  with  Impulse- 
2440  BSUBN=-ASUBN 

2450  BSUBP=ASUBP-2*RAT I 0*ASUBN* IMPSUBP 

2440  BSUBC2ASUBC-2*RATI 0*ASUBN* I MPSUBC 

2470  REM  GOTO  2710  s  REM  SKIP  GAUSSIttM  SPREAD 

2480  BOUNCES2" Bounce  with  Impulse  and  Gaussian  Spread- 

2490  IF  BOUNCE=0  THEN  3040 

2700  GOTO  2714  sREM  HIGH  DETUNING 

2702  SI GMAC2 . 2*SQR< BSUBN) 

2704  S I GMAP2 . 455*SQR <  BSUBN > 

2704  SI GMANI2 . 05*BSUBN+ .0000 1 2*BSUBNA3 
2708  SIGrSV^2SQR<SIGrS^PA2+SIGrW4IA2) 

2710  GOTO  2730 

2714  SI GMAC=.049*SQR< BSUBN) 

2718  S I GMAP2. 149*SQR< BSUBN) 

2720  SIGMANI2. 0 1 3*BSUBN+ .000001 4*BSUBNA3 
2722  S I GMAN*SQR <  S I GMAPA  2+ S I GMAN I A  2  > 

2730  TEMP2SQR<  2/PI )  ' 

2740  FOR  1=1  TO  3 
2750  PROX-RND 
2740  SPX-0 

2770  DSPX2.l*SGN<PR0X-.5> 

2780  DPF I R2 . 5-PROX 
2790  DPOLD-DPFIR 

2800  IF  N0T<SPX20>  THEN  DPOLD-DP 
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2810  SPX=SPX+DSPX 

2820  PROSP= . 5+TEMP* <  SPX/2-SPXA  3/12+  SPXA  5/80-SPXA  7/672+ SPXA 9/ 
69 12) 

2830  DP=PROSP-PROX 

2840  IF  DP*SGN<DPFIR> > .01  THEN  2800 
2850  IF  ABS(DP) < .001  THEN  2910 
2860  DSPX=-DSPX*DP/<  DP-DPOLD) 

2870  SPX=SPX+ DSPX 

2880  PROSP= . 5+TEMP*  <  SPX/2-SPX*  3/ 1 2+  SPXA  5/80 -SPXA 7/672+ SPXA  9/ 
6912) 

2890  DP=PROSP-PROX 

2900  GOTO  2850  s  REM  NEXT  ATTEMPT 

2910  IF  1=1  THEN  DGAUN=SPX 

2920  IF  1=2  THEN  DGAUP=SPX 

2930  IF  1=3  THEN  DGA(JC=SPX 

2940  NEXT  I 

2950  DN=S I GMAN*DGAUN 

2960  DP=SIGMAP*DGAUP 

2970  DC=SI GMAC*DGAUC 

2980  A*="SIGMAN=* 

2990  B*="VNORM= 

3000  PRINT  USING  F0RM2S ;A* ; SI GMAN; * DN=" ;DN; "DP=" ;DP; *DC=* ; 

DC ; B* ; BSUBN 

3010  IF  NOT  PRINTER=1  THEN  3030  :  REM  NO  PRINTER 

3020  LPRINT  USING  F0RM2*;A^;SIGMAN; *DN=" ;DN; "DP=B ;DP; BDC=" ; 

DC} 8*; BSUBN 
3030  BSUBN=BSUBN+  DN 
3040  BSUBP=BSUBP* DP 
3050  BSUBC=BSUBC+DC 
3060  RETURN 
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45  Degree  Incidence  to  Surface  Upward 
Bounce  with  Impulse  and  Gaussian  Spre. 
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